re Reta 
Balers aD 


Gx wpnis 
UNIOERSTTATIS 
AUBERTAEASTS 


Digitized by the Internet Archive 
in 2022 with funding from 
University of Alberta Library 


https://archive.org/details/Dingwell1984 


THE UNIVERSITY OF ALBERTA 


Investigations of the role of fluorine in silicate melts : 
implications for igneous petrogenesis. 
by 


Donald Bruce Dingwell 


A THESIS 
SUBMITTED TO THE FACULTY OF GRADUATE STUDIES AND RESEARCH 
IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE 


OF Doctor of Philosophy 


Department of Geology 


EDMONTON, ALBERTA 


Fall 1984 


ua); " . = pe coe Hey 

etl 

A Ma Mg 

a tus FY) 

en 
f o 
. 
: . yj 7 f] : 7 : , 
P wre 
i I ; 
arasigas 50 tans aur | | ti 


HORASAGR CMs 24 toute BRAUGAAD "FTO erswoKs ‘ 


bas (beer overt yous sere ehidiar benptesebay eft" 
foreseet bh6 esrbuse sipuberd fo y2iuont sit 02 
efor sai io anoitegitaievnl belirins ztasd? s ,@pns! tt 
audeipt tot anotis3i Lait’ + etion sisoilie “ snks 
fetiveq ai. f[! -eweinke soua@ ‘Bisnot ya bes? imdue., 2 
Lo ‘sotto, to -seapee od 02 ademas ivpe: ett * se i 7 


; ‘ey, 4 
ype f 
a , 
e ™/ . 
; i LP 
| | ee ey ; 
ae aes ch iediete le ie! an dN i ( 
i ' | ; 
‘ i) 
- . 
toe dv yegue” 7 7 
oe ore : oy | 


dis i ay 
b Aad 


Dedication 


This thesis is dedicated to my parents, Phyllis and Bruce. 


iv 


1 we ee _ a 


‘ "4 
¢ 1)" | al @ 


ABSTRACT 

The results of five separate studies on the properties of 
fluorine- and water-bearing silicate melts are presented in 
this thesis. Each study represents an attempt to evaluate 
the effects of fluorine on polymerized silicate melts by 
studying melt properties. The melt properties studied are 
viscosity, fluorine diffusivity, melt-vapor partitioning and 
water solubility. 

The effect of fluorine on melt viscosities of five 
compositions in the system Na,0-Al,0;-SiO, was investigated 
at one atmosphere and 1000-1600°C by concentric-cylinder 
viscometry. The compositions chosen were albite, jadeite and 
nepheline on the join NaAlO2-Si02z and two others off the 
join at 75 mole percent SiOz, one peralkaline and one 
peraluminous. Fluorine reduces the viscosities and 
activation energies of all melts investigated. The 
viscosity-reducing power of fluorine increases with the SiO, 
content of melts on the join NaAlO2-Si0O2 and is a maximum at 
Na/Al(molar)=1 for melts containing 75 mole percent SiOz. 
Fluorine and water have similar effects on aluminosilicate 
melt viscosities, probably due to depolymerization of these 
melts by replacement of Si-O-(Si,Al) bridges with Si-OH and 
Si-F bonds, respectively. The viscous flow of phonolites, 
trachytes and rhyolites will be strongly affected by 
fluorine. 

The chemical diffusion of fluorine in jadeite melt was 


investigated from 10 to 15 kbars and 1200 to 1400°C using 
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diffusion couples of jadeite melt and fluorine-bearing 
jadeite melt (6.3wt.% F). The diffusion profile data 
indicate that the diffusion process is 
concentration-independent, binary, F-O interdiffusion. The 
F-O interdiffusion coefficient ranges from 1.3 x 10-7 to 7.1 
x 10°’ cm?/sec. The Arrhenius activation energy of diffusion 
ranges from 36 to 39 kcal/mole as compared with 19 kcal/mole 
for fluorine tracer diffusion in a lime-aluminosilicate 
melt. The diffusivities of various cations are significantly 
increased by the addition of fluorine or water to a silicate 
melt. | 

The chemical diffusion of fluorine in five melts in the 
system Na2,0-Al203;-Si0O2 was investigated at 1200-1400°C and 1 
atm. Spheres of the same melts that were used for viscometry 
were suspended from Pt loops in a flow of oxygen gas, at 
high temperature. The resulting volatilization of fluorine 
was investigated by microprobe scans for Na, Al, Si and F, 
across quenched, sectioned spheres. Profiles of fluorine 
concentration can be reduced to yield a bulk diffusion 
coefficient for F-O interdiffusion. Fluorine diffusivity 
decreases along the join NaAl0O2,-SiOz from nepheline to 
albite. At 75 mole % SiO-,, fluorine diffusion is slower in 
albite than in peraluminous and peralkaline melts. Fluorine 
aiffusion in albite melt is significantly slower than water 
diffusion in obsidian melt. 

The solubilities of water in six melts in the system 


K,0-Na,0-Al,03-Si0O2 were determined at 930-1630 bars and 
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800°C. The solubilities were determined by micromanometric 
measurement of H20 in quenched melts. The solubilities at 
970 bars for the granitic and phonolitic minimum melts are 
2,00 20al0jwtthand 5S-0ivt0.14\ wt%;erespectively.nBoth 
peralkaline and peraluminous granitic melts have higher H,0 
solubilities than the 1 kbar P(H2O) minimum melt. Fluorine 
decreases the solubility of water in granitic melts. 

The partitioning of Na, K, Al, and Si between 
fluorine-bearing haplogranitic melts and a coexisting fluid 
phase was investigated at 1 kbar and 800°C. Fluorine 
Significantly increases the solubility of Na, K, Al and Si 
in the fluid phase. The largest increase is for aluminum 
suggesting that fluorine is coordinated by aluminum in one 
or more aluminofluoride or alkali aluminofluoride complexes 
under these conditions. The results provide experimental 
confirmation of the mobility of aluminum observed in many 


fluorine-rich magmato-hydrothermal systems. 
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I. Introduction 


A. Opening comments 

An accurate knowledge of the physical and chemical 
properties of silicate melts is essential for a discussion 
of the processes involved in igneous petrogenesis. This 
thesis contains five investigations into silicate melt 
properties which have direct consequences for the 
petrogenesis of many igneous rocks. Fluorine-rich igneous 
rocks occur in a variety of settings, both intrusive and 
extrusive, and the chemistry of fluorine-rich magmas can be 
wet or dry, peralkaline or peraluminous (Bailey, 1977). 
Typical and extreme ranges of fluorine contents of granitic 
and related rocks are shown in Figure 1. It is well known 
that fluorine significantly modifies several of the physical 
properties of synthetic aluminosilicate melts (see below) 
and, therefore, it is expected that fluorine-rich geological 
melts may behave differently than fluorine-free melts of 
equivalent composition. This thesis comprises a series of 
studies which investigate the role of fluorine in silicate 
melts of geologic interest by examining some of the chemical 


and physical properties of fluorine-bearing melts. 


B. Background 


Previous work on the role of fluorine in silicate melts 


can be grouped into two categories. Firstly, there has been 


a long discussion in the ceramics literature on the role of 
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Figure 1. Common and extreme ranges of fluorine contents of 
granites and related rocks. 
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fluorine in silicate melts and glasses at one atmosphere. 
Studies have been conducted on viscosity (Schwerin, 1934; 
Owens-Illinois Glass Company, 1944; Kozakevitch, 1954; 
Bills, 1963; Hirayama and Camp, 1969; Rau et a/., 1977), 
electrical conductivity (Shinozaki ef/ai.)) 1997) eeintrared 
and Raman spectroscopy (Kumar et ai., 1965; TakuSagawa, 
1980; Dumas et a].; 1982), volatilization (Kumar et al]., 
1961; Kogarko et a/j., 1968; Al-Dulaimy, 1978), melting 
curves (Kogarko and Krigman, 1973) and opalescence 
(Rothwell, 1956). With the exception of some of the 
viscosity work, all these studies have been conducted on 
Gepolymerized, Al-poor melts. The major conclusion to be 
drawn from these ceramic studies is that fluorine is 
incorporated in dry, aluminum-free melts by forming Si-F 
bonds, replacing the oxygen in Si-O-Si bridges, and 
resulting in the formation of Si0O;F units (Rabinovich, 
1967). Rabinovich (1983) has discussed the consequences of 
structural reorganization during slow cooling of 
fluorine-bearing silicate melts (i.e. opalescence) and the 
increasing volatilization of fluorine as alkali fluorides 
with decreasing temperature. Studies which suggest the 
predominance of alkali-fluoride bonds in Silicate melts 
(Kogarko et al., 1968; Kogarko and Krigman, 1973) have 
probably undergone the type of structural reorganization 
described by Rabinovich (1983) or have been conducted on 


extremely alkaline, depolymerized melts. 
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The second body of data on fluorine in silicate melts 
comes from extensive phase equilibria work on wet, 
fluorine-bearing granitic and albite melts (Wyllie and 
Tuttle, 1961; Koster van Groos and Wyllie, 1968; Anfilogov 
et al]l., 1973; Glyuk and Anfilogov, 1973a,b; Danckwerth, 
1980; Manning, 1981). Despite the variety of synthetic and 
natural compositions used, almost all studies report the 
replacement of feldspar by quartz on the liquidus with 
sufficient addition of fluorine. This relative increase in 
quartz stability implies a decrease in the activity of 
NaAlO,, relative to $102, in the melt. Decreased NaAlO, 
activity is not consistent with the solution of fluorine in 
these melts by the formation of Si03F complexes. Instead, 
Manning et a]. (1980) have proposed the formation of Al-F 


complexes with Al in octahedral coordination. 


C. Objectives 

This thesis is an attempt to evaluate the behavior of 
fluorine and fluorine-bearing melts in simplified chemical 
systems. The results of these studies are discussed in terms 
of their implications for the structural role of fluorine in 
Silicate melts and the petrogenesis of fluorine-bearing 
igneous rocks. The thesis is written in the form of five 
separate papers (chapters 2 to 6) and therefore each chapter 


has its own introduction and conclusions. The conclusions of 


all studies are summarized in chapter 7. 
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Viscometry 


Chapter 2 is an investigation of the effect of fluorine 
on melt viscosities in the system Na,0-Al,03;-Si02. Viscosity 
iS a property which has been shown to vary greatly with melt 
composition in simple systems (e.g. Riebling, 1966) and 
which has been related to melt structure (Mysen et al., 
1980). The viscosity of igneous magmas is a critical 
parameter controlling their physical behavior during 
coalescence, ascent and emplacement/eruption (Harris et al., 
1970). The primary objective of this work was to quantify 
the viscosity-reducing effect of fluorine on relatively 
pewmenized melts in the system Na20-Al203;-Si0O,2. The results 
indicate that the viscosity-reducing effect of fluorine is 
comparable to that of water. Fluorine and water probably 
remain stably dissolved in silicate melts by replacement of 
PaineSs-O40Sij}Al) wbridges with (Si,Al)-F and (S1,A1)-0H, 
respectively. The viscosity-reducing effect of fluorine is 
greatest in SiO2-rich melts at low temperatures (600-800°C) 
where viscosity decreases several orders of magnitude with 


thenadartroncofe teweseftluorine: 


Diffusion 

Chapters 3 and 4 present data on the diffusivity of 
fluorine in melts of geologic interest. Diffusivities in 
silicate melts yield information on the structural nature of 
ionic transport in melts which can be usefully compared with 


viscosity and other melt properties to obtain a clearer 
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understanding of melt structure. The geologic applications 
of fluorine diffusivity concern bulk diffusion rates in 
relatively dry, fluorine-rich igneous magmas. 

The chemical diffusion of fluorine has been studied at 
1 atmosphere in five melts in the system Na,0-Al,0;-Si0O, and 
at 10 to 15 kilobars in jadeite melt. The diffusion is a 
binary F-O exchange which is independent of concentration 
but strongly dependent on melt composition. Both fluorine 
and water significantly affect the diffusivities of other 
species in silicate melts and thus the bulk diffusivities of 
magmas are enhanced by their presence. There are several 
examples of relatively dry, fluorine-rich geologic melts 
whose kinetic behavior may be controlled their fluorine 


contents: 


Water solubility 

Few properties of silicate melts are more influential 
in igneous petrogenesis than the solubility of water. Yet 
the data available for several natural melts contain 
Significant discrepancies. Chapter 5 presents the results of 
water solubility determinations in haplogranitic melts using 
a new technique, vacuum fusion manometry. The solubility of 
water was measured in 1 kilobar thermal minimum, peralkaline 
and peraluminous granitic melts, in a phonolite melt and in 
fluorine-bearing granitic melts. Water solubility varies 
strongly with melt composition. The solubility of water in 


peralkaline and peraluminous melts is greater than in the 1 
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kilobar minimum granitic melt and the phonolite melt 
dissolves twice as much water as granitic melts. The 
addition of fluorine to the 1 kilobar minimum melt decreases 


the solubility of water slightly. 


Melt-fluid partitioning 

The partitioning of elements between a melt and a 
coexisting fluid phase is a critical parameter in the 
evolution of magmato-hydrothermal systems. 
Chlorine-dominated systems are common but fluorine-rich 
systems do exist. Chapter 6 reports oar eerie data for 
Na, K, Al and Si in the system haplogranite-HF-H,0. Fluorine 
enhances the solubility of Na, K, Al and Si despite the high 
melt/fluid partition coefficient of fluorine itself. 
Normalized fluid compositions in equilibrium with 
fluorine-bearing granitic melts are much more Al-rich than 
comparable chlorine-bearing and boron-bearing fluids. 
Fluorine does not supress Al mobility in 
magmato-hydrothermal systems as chlorine does and therefore 
fluorine-rich systems can be expected to, and often do, 


exhibit evidence of extensive Al mobility. 


Summary 


Chapter 7 is a brief summary of the geological and melt 


Structural implications of the contents of chapters 2 to 6. 
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II. Viscometry 


A. Introduction 

It is well known in the glass and ceramics literature 
that the incorporation of fluorine into silicate melts has 
an important effect on their properties (Weyl, 1950; Eitel, 
1965). More specifically, fluorine, as CaF, has been used 
in ceramics for several hundred years as a flux to promote 
fining or bubble removal from glass melts and more recently 
as a flux to aid the kinetics of slag-iron separation in 
steelmaking. The ceramics and glass literature contains a 
large body of data on the effect of fluorine additions to 
various slag compositions (Schwerin, 1934; Owens-Illinois 
Glass Company, General Research Laboratory, 1944; 
Kozakevitch, 1954; Bills, 1963; Hirayama and Camp, 1969). In 
the geological literature it has been shown that fluorine 
plays an important role in melt structure and properties 
(Wyllie and Tuttle, 1961; Manning, 1981). 

For these reasons the effect of fluorine on the 
viscosity of melts in the system Naz0-Al203;-Si0O2 was 
investigated. Temperature-viscosity relationships of 
fluorine-bearing melts of albite, jadeite, nepheline, one 
peralkaline and one peraluminous composition have been 
measured between 1000-1600°C. The results were then compared 
with available data from fluorine-free melts in 


Na,0-Al,03;-Si0O2, with studies on fluorine-bearing slags, and 
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with data for hydrous rhyolites. From these comparisons it 
PSeriterred ‘that (1) reductions invviscosity of polymerized 
melts occur by substitution of fluorine for bridging oxygens 
and (2) water and fluorine have similar effects on the 


viscosities of polymerized melts. 


B. Experimental method 

The compositions studied were chosen to evaluate the 
influence of varying alkali/aluminum ratio and silica 
content on the viscosity of fluorine-bearing melts. The five 
base compositions are albite, jadeite, nepheline and two 
others with 75 mole percent SiOz, one peralkaline and one 
peraluminous (Fig. 2). The system Na20-Al203;-Si0O2 was chosen 
because viscosities in the base system Na;0-Al20;-Si02 have 
been studied by Riebling (1966) and the structure and 
physical properties of glasses and melts in this system have 
been investigated (Hunnold and Bruckner, 1979; Taylor and 
Brown, 1979;"Navrotsky et al., 1982: Seifert ef a/., 1982). 
The normative compositions of many fluorine-rich igneous 
melts containing greater than 50 mole percent feldspars + 
feldspathoids are well represented by this system. 

Fluorine waS included in the melt compositions by 
substituting 2AlF; for some of the Al20; of the base 
composition. Thus the substitution is essentially 2 moles of 
fluorine for one mole of oxygen, denoted by the exchange 
Operator F,0.,. The starting compositions were synthesized 


from reagent grade sodium carbonate, alumina, aluminum 
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Figure 2. Melt compositions projected into Na20-Al203;-Si02. 
One atm liquidus phase fields; cristobalite(cr), 
tridymite(tr), albite(ab), sodium disilicate(Na-di), 
sodium metasilicate(Na-mt), nepheline(ne), 
carnegieite(cn), corundum(co) and mullite(mu). The 
diagram is in weight percent. 
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1 atm 
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~peraluminous’ 


* peralkaline’ 
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fluoride and purified quartz sand. Carbonate + oxides + 
fluoride, equivalent to a decarbonated weight of 600 grams, 
were mixed thoroughly for 12 hours and then fused in a 10.5 
cm long by 5.5 cm diameter platinum crucible for 6 to 10 
hours at 1600°C to ensure homogeneity and escape of air 
bubbles. 

Viscosities were determined using a concentric cylinder 
viscometer which uses a Pt90Rh10 inner cylinder 5 cm long 
and 1.2 cm in diameter with conical ends. The inner cylinder 
is rotated and the resultant torque is measured and 
converted into a millivolt signal. The apparatus was 
calibrated with NBS standard lead-silica glass SRM #711 for 
which the viscosity-temperature relationship is well known. 
The viscosities are accurate to +5% with a precision of +1% 
and temperatures have uncertainties of +1°C. Measurements 
were made at 50°C intervals over the temperature range 
1000-1600°C. Viscosities were independent of shear rate for 
all compositions (e.g. Fig. 3) indicating Newtonian 
behavior. 

Volatilization of fluorine was significant during melt 
synthesis but not during viscosity measurement. The 
viscosity measurements were obtained at successively lower 
temperatures for each composition and then the highest 
temperature measurement was repeated. Due to the powerful 
viscosity-reducing effect of fluorine in these melts, 
returning to the high temperature data point for each 


composition is an excellent check for significant fluorine 
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Figure 3. Newtonian behavior of fluorine-bearing (5.8 wt.%) 


albite melt (data at 1500°C). 
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volatilization during measurements. High temperature 
readings at the start and finish of each set of measurements 
were within error, indicating that fluorine volatilization 
during the measurements was insignificant. According to the 
glass compositions analyzed after viscometry, there was no 
Significant sodium loss. 

Na, Al and Si contents of quenched melts (glasses) were 
determined, after completion of the viscosity measurements, 
by energy dispersive analysis using an ARL-SEMQ microprobe 
fitted with an EEDS-ORTEC energy dispersive system. 
Operating conditions were 15 kV, 4nA sample current and 240 
second count times. The beam was rastered over a 20x20 
micrometer area, a technique which proved adequate for 
avoiding volatilization of Na or F during analysis. The 
homogeneity of glasses was confirmed for Na, Al and Si by 
analysing six spots on each glass. 

Fluorine contents were determined by neutron activation 
analysis. Twelve replicates of each glass were determined 
against a curve for reagent-grade CaF2. The technique was 
verified with standard opal glass SRM #91 (5.72 wt.“F). 
Fluorine contents were combined with the raw spectra for Na, 
Al and Si as input for EDATA2 (Smith and Gold, 1979), 


allowing full ZAF reduction of the Na, Al and Si data (Table 
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Table 1. Analyzed melt compositions* 


Na Al on F O Total 
ilbite seal 0.15 Site 4 Se! 44.5 99.30 
jadeite 1B leaping se as) ZO a02 Ons a 8) 99.94 
iepheline ESO 7 Los 19.40 Dae Alt 99.24 
yeralkaline 10.88 42 34.09 atl 45.9 ea ehe, 
yeraluminous G.14 10.87 Sm, 5S 5.4 AD 99.49 


‘Na, Al and Si determined by electron microprobe; F determined by neutron 
ictivation analysis; O by stoichiometry. 

irrors for microprobe data expressed as percent of the amount present at 3 
mandard-deviations: »Na (4.34), Al «(2<04) ,f SiC! .OZ):. 

irrors in fluorine determinations are +0.lwt% F at 1 standard deviation. 
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C. Results 

Fifty viscosity measurements were made on five 
compositions. The results are presented in Table 2 and in 
Figure 4. All five compositions show a log-linear dependence 
of viscosity on reciprocal temperature. Such behavior is 


described mathematically by the following Arrhenius equation 


10g, 9 = 109) oN + En/2 .303RT 


where n is the viscosity at temperature T (in K), Ho is a 
constant, R is the gas constant and En is termed the 
activation energy of viscous flow. The method of least 
Squares was used to fit the data to straight lines (Fig. 4) 
from which the Arrhenian parameters were determined. Above 
1400°C, viscosities of the fluorine-bearing melts decrease 
in the order peraluminous > albite > jadeite > peralkaline > 
nepheline. Activation energies of viscous flow of 
fluorine-bearing melts decrease with increasing $102 along 
the join NaAlO,-SiO, from nepheline (58.0 kcal/mole) to 
albite (45.5 kcal/mole), and increase from the peralkaline 
melt (40.3 kcal/mole) to the peraluminous meit (68.7 
kcal/mole) at 75 mole percent S$iO,. Different activation 
energies for the various melts result in two 
viscosity-temperature curve crossovers within the range of 
experimental measurements (Fig. 4). As a result, below 
1350°C fluorine-bearing jadeite is more viscous than 


fluorine-bearing albite, and below 1400°C fluorine-bearing 
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Figure 4. Measured viscosities of fluorine-bearing melts 
(this study) and base melts (Riebling, 1966). Numbers on 
curves refer to the wt% of fluorine. Inset: Viscosities 
of vitreous and fluorine-bearing SiO, (Rau et al., 
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Table 2. Viscosity results and computed Arrhenian parameters 


ea Sth SSD 
a ea ec a a ee a Se 


Temperature Log}9 Viscosity (poise) 
(degrees C) albite jadeite nepheline peralkaline peraluminous 
1600 Ld 6 2.45 - - S12 
1550 2.93 2.08 1.37 2.04 Sivas 
1500 3.10 2.88 2.10 Zea” 3.50 
1450 5.28 3.08 fel 2632 3.02 
1400 3.46 3.29 2.03 2.44 4.09 
1350 3.00 Bool PBING Pan S) 4.38 
1300 3.84 eho cael) Ziel 4.68 
1250 4.05 3797 S226 2.97 4.96 
1200 4.25 ae - 346 Fe28 
1150 4.50 G62” - S2e7 - 
1100 LAST 7 - 3.60 = 
1050 4,96 ~ = 3.86 = 
1000 2 - ~ 4,14 = 
Arrhenian parameters 
ae 45.5 55.0 58.0 40.3 68.7 
logio99 -2.50 -3.92 -5.05 -2.81 -4,89 
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nepheline is more viscous than fluorine-bearing peralkaline 


melt (Fig. 4). 


D. Discussion 


Comparison with previous work. 

Discussion of these data is facilitated by comparison 
with viscosities of equivalent fluorine-free melts in 
Na20-Al203;-Si0O2. In Figure 4 we have included data from 
Riebling (1966) for the five base compositions and the data 
of Rau et al. (1977) for vitreous Silica and 
fluorine-bearing Silica. Figure 4 shows that the addition of 
fluorine strongly reduces the viscosity and the activation 
energy of all melts studied. These reductions 
(denoted dn and dEn, respectively) are different in 
magnitude for each of the base compositions. The lowering of 
viscosities and activation energies of melts on the join 
NaAlO.,-SiO, is a positive function of S102 content. For 
example, at 1400°C the viscosity reduction with the addition 
of 1 wt% fluorine (by substitution of fluorine for oxygen) 
£otSiOsiis 0745el0g,$uunits: | Undersequivalent conditions, 
the viscosities of albite, jadeite and nepheline melts 
decrease 0.26, 0.19 and 0.16 logio units/wt.% F, 
respectively. All activation energies for fluorine-bearing 
melts are reduced from the values for their fluorine-free 


counterparts: albite (85.3 kcal/mole), jadeite (86.8 
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kcal/mole) and nepheline (74.2 kcal/mole) (Riebling, 1966). 
Interestingly, the decrease in activation energy of these 
three melts upon addition of fluorine is a linear function 
of mole fraction of Si02/(Si0,+NaAl0,) over the range of 0.5 
to 0.75 (Fig. 5). This results in decreasing activation 
energies with increasing SiO, content from fluorine-bearing 
nepheline to fluorine-bearing albite. Extrapolation of the 
finé in Figure 5 to X(Si0O.)=1 is not possible due to the 
large uncertainty in the activation energy of silica (see 
Hofmaier and Urbain, 1968) and the possibility that 1 and 6 
wt% additions of fluorine may not have an equivalent effect 
on the activation energy of S$i0,. However, the activation 
energy decrease is larger for the addition of fluorine to 
Silica than for the addition of fluorine to any Al-bearing 
melts. 

The three melts studied at 75 mole percent SiO2 are: 
peralkaline, albite and peraluminous. At 1400°C, the albite 
melt undergoes the largest decrease in viscosity (0.26 log;o 
units/wt.%F), whereas the peralkaline and peraluminous melts 
show smaller decreases of 0.16 and 0.14 logio unitS/wt.4%F, 
respectively. These three melts all experience decreases in 
activation energy with fluorine addition, but the 
peralkaline and peraluminous melts drop only 5 and 16 
kcal/mole, respectively, compared to a 40 kcal/mole drop for 
pee eiAs a result the viscosities and activation energies 


increase in the order peraluminous > albite > peralkaline. 
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Figure 5. Decreases in activation energy of viscous flow 
(den) aS a function of mole fraction .of /Si0, (Ret 
al.=Rau et a]., 1977; H=Hofmaier and Urbain, 1968). 
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In summary, in order to account for the role of 
fluorine in melts in the system Na20-Al203;-Si02 any model 
must explain the following general features: 1) the decrease 
of all viscosities and activation energies with fluorine 
addition; 2) the positive dependence of dEn and dn (see 
above) on SiOz content for melts on the join NaAlO.-SiO., and 
the resulting negative dependence of En on SiO, for 
fluorine-bearing melts on this join; and 3) the positive 
dependence of En and on on Al/Na at constant SiO2, that 
contrasts with the presence of viscosity maxima at or near 
Al/Na (molar)=1 for fluorine-free melts (Riebling, 1966). 

There have been several studies conducted on the 
effects of fluorine on viscosities of Al-poor slags 
(Schwerin, 1934; Owens-Illinois Glass Company, General 
Research Laboratory, 1944; Kozakevitch, 1954; Bills, 1963; 
Hirayama and Camp, 1969). The compositions of slags vary 
widely, especially in silica content (25 to 75 wt%) but, due 
to insufficient aluminum, none of the melts represents 
highly polymerized liquids such as those of the join 
NaAlO.-SiO,. Representative compositions and results from 
these studies are presented in Figure 6 with the results of 
the present study, thus providing a comparison of the effect 
of fluorine on polymerized and depolymerized melt 
viscosities. It is clear in Figure 6 that fluorine decreases 
the viscosities of all melts studied. However, the decreases 
in activation energy with fluorine addition are smaller for 


the depolymerized slags than for the polymerized melts in 
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Figure 6. Effect of fluorine on melt viscosity for several 
Slag compositions. See text for references. Numbers on 
curves refer to wt.% fluorine added. 
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Na 20-Al 203-S103 ; 


Structural implications 

The three-dimensional structure of melts on the join 
NaAlO2-Si0O2 has been discussed by several investigators 
(Riebling, 1966; Hunnold and Bruckner, 1979; Taylor and 
Baownyi. 19795 iNavrotsky et.al 19° 1982eeSeifertneteal .-ni982). 
A continuous three-dimensional network of (Al,Si)O, 
tetrahedra is stabilized by the inclusion of one network 
modifying Na atom per tetrahedrally coordinated Al atom. In 
such a structure all oxygen is bonded to two 
tetrahedrally-coordinated, network-forming cations to form 
oxygen bridges. In contrast, melts in Na20-Al,03;-Si0, whose 
compositions lie off the join NaAlO2-Si0.2, whether they are 
peralkaline or peraluminous, cannot maintain the structure 
of an uninterrupted three-dimensional network. These melts 
probably contain a discontinuous three-dimensional network 
perturbed by the presence of excess Na or Al atoms, 
coordinated by nonbridging oxygens (NBO) (Mysen et al., 
1980). The degree of polymerization of aluminosilicate melts 
may be represented by the ratio of NBO/T (where T represents 
a tetrahedrally-coordinated cation ; Mysen et al]., 1982). 
Analyses of the melts along the join NaAlO2-SiOz yield base 
compositions (projected from F,0.,(molar)) with NBO/T values 
close to zero (albite = 0.013; jadeite = 0.008; nepheline = 


0.032). The peralkaline base composition yields NBO/T4= 
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0.193. The peraluminous composition is difficult to deal 
with using the calculation procedure of Mysen et al. (1982) 
because of the uncertainty in assigning some or ail of the 
aluminum to the value of T. If we assume as a limiting case 
that each tetrahedrally coordinated Al atom requires a Na 
atom for stabilization, then the maximum NBO/T for the 
peraluminous composition is 0.292. Viscosities in 
Na20-Al20;-Si0O2 (Riebling, 1966) appear to reflect the 
Structural characteristics of these melts quite well. The 
presence of viscosity maxima at or near Na/Al (molar) = 1 
along joins of constant silica content in Na2,0-Al,0;-Si0>; 
(Riebling, 1966) imply that at the 1:1 composition melts are 
completely polymerized (i.e. all Al is in tetrahedral 
coordination; Hunnold and Bruckner, 1979). 

The viscous flow of pure silica involves an activation 
energy similar to the Si-O bond strength of Si-O-Si bridging 
bonds and thus viscous flow in SiOz is thought to proceed by 
the breakage of such bridging bonds (Bockris and Reddy, 
1970). Melts along the join NaAl02-SiOz also require the 
breakage of bridging oxygen bonds but the average T-O-T bond 
strength is lowered from that of pure SiOz. Peralkaline and 
peraluminous melts in Na2,0-Al203;-Si0Oz probably owe their 
lower viscosities to the smaller average size of flow units. 
Their lower activation energies may be due to a decrease in 
the number of bond breakages required for viscous flow. 

Fluorine substitutes for oxygen in silicate melts 


(Rabinovich, 1983). If oxygen occurs in melts on the join 
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NaAlO2-SiO2 entirely as Si-O-(Si,Al) bridges, then fluorine 
must break oxygen bridges. It is clear from the viscosity 
data for fluorine-bearing silica (Rau et al., 1977) that the 
exchange of Si-O-Si bridges for Si-F bonds occurs with a 
drastic change of melt viscosity (Fig. 4, inset). However, 
for melts along the join NaAlO2-SiO,z the question is whether 
Al-O or Si-O bonds, or both, are the primary target for 
fluorine substitution. Fluorine-free melts have decreasing 
viscosities in the order albite > jadeite > nepheline, but 
Similar activation energies (Riebling , 1966). The magnitude 
of the effect of fluorine on both activation energies and 
viscosities of melts on the join NaAlO2-SiO2z iS a positive 
function of SiOz content. Figure 5 therefore supports the 
proposal that the decreasing efficiency of activation energy 
reduction with decreasing SiOz content is a result of the 
decrease in the proportion of Si-O-Si bridges in these 
melts. 

The three melts at 75 mole percent SiOz vary in Al/Na 
and their viscosities and activation energies increase with 
Al/Na. Apparently the addition of fluorine to melts along 
this join eliminates the viscosity maximum at or near 
Al/Na=1 for fluorine-free melts. Fluorine in the peralkaline 
and peraluminous melts may be substituting for non-bridging 
or bridging oxygens. 

Breakage of Si-O-Si bridges can also explain the effect 
of fluorine on the viscosities of depolymerized melts such 


as those included in Figure 6. Melts shown in Figure 6 have 
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NBO/T ranging to 4.7 and thus are more depolymerized than 
the melts investigated in this study. Despite their 
calculated values of NBO/T, the SiO2-poor melts probably 
contain some bridging oxygens (Gotz et al., 1976; Smart anda 
Glasser, 1978). Reduction of viscosities by fluorine 
addition to these melts may result from substitution of 


fluorine for bridging or nonbridging oxygens. 


Comparison with water 

Water has a strong influence on the viscosity of 
Silicate melts (Shaw, 1963; Friedman et a/]., 1963; Burnham, 
1964: Scarfe, 1973). The similar size and charge of F~ and 
OH” ions encourage comparison of their effects. We would 
prefer to compare the influences of water and fluorine on 
identical melt compositions; however, the data do not exist. 
Figure 7 contains viscosity-temperature curves for dry and 
wet rhyolites (Shaw, 1963; Friedman et a]J., 1963; Carron, 
1969), a wet pegmatite (Burnham, 1964), albite (Riebling, 
1966) and fluorine-bearing albite (this study). Activation 
energies decrease similarly with addition of fluorine to 
albite and addition of water to rhyolite. This implies 
Similar roles for fluorine and water in the depolymerization 
and resulting reduction of viscosity in these polymerized 
melts (Fig. 7). The magnitude of viscosity reduction is a 
much more difficult comparison. If we select a dry rhyolite 


curve to coincide with the albite curve then the reduction 
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Figure 7. Comparison 
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of viscosity (normalized on a weight percent or mole percent 
basis of added water or fluorine) is slightly greater for 
the addition of water to rhyolite than for the addition of 
fluorine to albite. However, recalling the positive 
relationship between SiO, content and the viscosity-reducing 
power tofmiluorines(PigeesS);sit isllrkely thatWcompanison lof 
two rhyolites rather than albite and rhyolite would remove 
this discrepancy (note that converting from weight to mole 
percent does not significantly alter comparisons between 
Fluorine.and water as their molecular weights are 19.0 and 
18.0 grams/mole, respectively). Therefore, the conclusion 
from comparison of water and fluorine is that one mole of 
fluorine is structurally analogous to one mole of water in 
depolymerizing highly polymerized melts. The equivalence of 
the viscosity reductions strongly implies that H20 and 
fluorine disrupt an equivalent number of oxygen bridges when 
mcosporatediim highly spolymerrzedcmelts.UThe eftecttot 
fluorine on hydrous melt viscosities cannot be discussed at 
this time because, with the exception of the Harding 
pegmatite (0.33 wt.% F, 8.8 wt.%H20), there are no data on 


the viscosities of fluorine-bearing hydrous melts. 


Geological applications 
Many volcanic rocks have anhydrous chemical 
compositions very close to the compositional plane 


KA1Si0,-NaAlSi0,-SiO,. In particular, "average" compositions 
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of phonolites, trachytes and rhyolites (Cox et a/]., 1979, 
p.401) have values of K20+Naz0+Al1203+Si02 near 90 weight 
percent. Thus the bulk physical properties of melts in the 
System K,0-Na20-Al203;-Si02, closely resemble the properties 
of these natural melts. Also, the molar effects of KO and 
Naz2O on melt viscosities are similar and K,0+Na-O (molar) is 
a single parameter in Shaw's (1972) scheme for calculation 
of natural melt viscosities. The molar equivalence of Na,O 
and K20 means that the system K20-Na20-Al,03;-SiO,2 may be 
effectively reduced to the simpler system Na,0-Al,0;-Si0O, 
that 1s used as a basis for this investigation. These 
arguments allow application of these results to natural 
melts. In the previous discussion, the importance of SiO, 
content in determining the ability of fluorine to reduce 
viscosities was discussed. The observation that the 
magnitude of the reduction in viscosity is positively 
correlated with $i0O2 is illustrated in Figure 8, where logio 
viscosity is plotted versus mole fraction SiO, for several 
compositions. Rather than plot a mixture of experimental and 
calculated viscosities in Figure 8, all viscosities were 
calculated using the method of Shaw (1972). Figure 8 
contains melt compositions expressed in terms of their mole 
fraction of $iO,. All natural and synthetic melts were cast 
into desilicated norms (composed of SiOz, NaAlO;z, CaAl20,, 
FeO, MgO, etc.) from which mole fractions of Sid, 
(Si0,/(Si0,+NaAlO2+KA10,+CaAl20,4,etc.)) were obtained. The 


synthetic phonolite, trachyte and rhyolite are anhydrous 
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Figure 8. Application of experimental results to synthetic 
The effect of fluorine content on 
melt viscosity for rhyolites, trachytes and phonolites 


and natural melts. 


at eruption temperatures near 1000°C. See text for 
explanation. 
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projections of the 1 kbar thermal minima in the system 
KA1S10,-NaA1Si0,4-SiO,-H20 (Tuttle and Bowen, 1958: Hamilton 
and MacKenzie, 1965) and the three natural compositions are 
maverages” +fromeCox et ali; y(1979)) Thetlines@labellede) 12) 
3, 4 and 5 represent viscosity decreases with weight percent 
additions of fluorine that are derived from linear 
interpolation of the data of this study. At an eruption 
temperature of 1000°C the viscosities of phonolites, 
trachytes and rhyolites containing 1 weight percent fluorine 
will be approximately 0.25, 0.5 and 1 orders of magnitude 
lower than the viscosities of equivalent fluorine-free 
melts. In comparison, steshaly peraluminous or peralkaline 
melts will experience a smaller influence of fluorine on 
their viscosities. 

The most dramatic effects of fluorine addition will be 
felt in relatively dry, fluorine-rich, high-silica 
rhyolites. Pre-eruption water contents of many silicic 
magmas may be lost during degassing immediately prior to 
eruption (Sparks, 1978). Due to a high melt/fluid partition 
coefficient (Hards, 1978), fluorine will be retained in the 
melt. Fluorine remaining in erupted lavas will maintain low 
viscosities. In certain cases fluorine may significantly 
influence eruptive style enhancing fluid flow and strongly 
reducing the pyroclastic component of a given eruption. 


Fluorine-rich topaz rhyolites exhibit just such features 


(Christiansen et a]., 1983). 
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Fluorine-rich intrusive rocks occur in a wide variety 
of settings. The rates of crystal settling and bubble ascent 
in such intrusions will be enhanced by high fluorine 
contents. The relatively dry magmas which form fluorine-rich 
alkalic intrusives will experience reduced viscosities as 
outlined above. The water-saturated emplacement of 
fluorine-rich melts such as Li-F leucogranites and 
associated pegmatites (Bailey, 1977; Manning, 1981) may be 
Significantly more fluid due to an additive effect of 
fluorine and water on the viscosities of such melts. At the 
greater depths and lower temperatures associated with 
water-Saturated emplacement of leucogranitic melts the 
distinction between low viscosity fluorine-rich, 
water-Saturated melts and coexisting fluids may be 
diminished. This distinction is the subject of considerable 
controversy regarding the magmatic or hydrothermal origin of 
certain late-stage fluorine- and water-rich cupolas and 


dykes (e.g. Kovalenko, 1973; Eadington and Nashar, 1978). 


E. Summary and conclusions 

It has been shown that the addition of fluorine to 
melts in the system Na2,0-Al,0;-Si02 strongly reduces 
viscosities and activation energies. Viscosities of 
fluorine-bearing aluminosilicate melts are Arrhenian 
functions of temperature. The viscosity-lowering effect of 
fluorine increases with increasing SiOz content along the 


join NaAlO,-SiO, and is a maximum at or near Na/Al (molar) 
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sie totemelts containinge?s pmoleo%eSi0s futheumost probable 
explanation for viscosity decrease is depolymerization of 
melts by replacement of Si-O-(Si,Al) bridges with Si-F 
bonds. When the effects of water and fluorine are compared, 
Similar reductions in viscosities and activation energies 
suggest Similar mechanisms of melt depolymerization. The 
most important implications of this work concern relatively 
dry, crystal-poor, magmas of phonolite, trachyte and 
rhyolite composition and their intrusive equivalents, both 
wet and dry. Due to the activation energy decreases, the 
viscosity-lowering effect of fluorine will be greater at one 
lower temperatures associated with water-saturated 
emplacement of igneous melts at depth. It is possible that 
the effects of fluorine and water might be additive, 
resulting in extremely large viscosity reductions for 


certain melts. 
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III. High pressure diffusion 


A. Introduction 

A knowledge of the transport properties of melts of 
geologic interest is required in order to model their 
behavior during petrogenetic processes. In particular, 
cationic and anionic diffusivities in silicate melts allow 
us to relate time, temperature, and physical scale in 
processes where diffusion is the rate-limiting step. 
Examples of such igneous processes include both melt-vapor 
and melt-crystal interactions (e.g. growth BE zoned 
minerals, vapor phase transport of dissolved metals during 
magma degassing, crustal assimilation) and intra-melt 
processes such as thermogravitational diffusion and 
double-diffusive convection. 

Fluorine, like water, has a considerable effect-on many 
properties of silicate melts including viscosities (chapter 
1), phase equilibria (Manning et a]., 1980), melt-vapor 
partitioning (Hards, 1978), and, as discussed below, 
component diffusivities in the melt. The occurrence of 
fluorine-rich amphiboles and micas in the lower crust and 
upper mantle (Smith et al., 1981; Valley et a/J., 1982) and 
the suggestion that many relatively dry, fluorine-rich, 
melts originate from these regions (Harris and Marriner, 
1980; Burt et al]., 1982) indicate the need for a better 


understanding of the role of fluorine in melts at high 
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pressures. Such considerations, along with the potential for 
insights into the structure of F- and H2O-bearing melts, 
prompted this study. Jadeite melt was chosen for this study 
because it has been used as a model for polymerized silicate 
melts in several studies (e.g. melt viscosity, Kushiro 
(1976); oxygen diffusivity, Shimizu and Kushiro (1984); 
cationic diffusivities, Kushiro (1983); Raman spectra, 
Sharma et a]. (1979), Seifert et a]. (1982); 


fFluorine-bearing melt viscosity, chapter 1). 


B. Experimental method 

The starting materials were (1) fluorine-bearing 
jadeite glass prepared from reagent-grade sodium carbonate, 
alumina, aluminum fluoride and purified quartz sand and (2) 
jadeite glass prepared from a gel, dehydrated at 800°C for 
two hours. These starting glasses were analyzed for Na, Al 
and Si by energy dispersive analysis using an ARL SEMQ 
microprobe fitted with an EEDS-ORTEC energy dispersive 
Spectrometer. Operating conditions were 15 kV accelerating 
voltage, 4 nA sample current and 240 sec counting time. The 
fluorine content of the fluorine-bearing starting glass was 
determined by neutron activation analysis as described in 
chapter 1. Analyses of the starting glasses are presented in 
Paure ss. 

The diffusion couple technique of Kushiro (1983) was 


used for this study. Glasses were ground in an agate mortar 
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Table 3. Analyses of starting glasses* 


A Sr Snr 
SE 


Element Jadeite Fluorine-bearing Stoichiometric 
Jadeite Jadeite 
Na d1ey5 PP PEC) 
Al 134.34 12275 13.34 
Si 27.92 eo de 27-79 
0 47.56 43.00 47 .54 
F - 6.3 - 
TAL 99.97 99.94 100 .00 


3, Al and Si determined by electron microprobe; F determined by 


Itron activation analysis; 0 by stoichiometry. 

rors for microprobe data expressed as percent of the amount present 
3 standard deviations: Na (4.3%), Al (2.0%), Si (1.0%). 

rors in fluorine analysis are + 0.1 wt% F at 1 standard deviation. 
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and packed into (5 mm diameter by 8-10 mm length) platinum 
capsules using a tight-fitting, stainless steel tool. The 
denser, fluorine-free powder was packed into the capsules 
first and occupied the lower end of the vertical diffusion 
couple in all experiments. Packed and crimped capsules were 
dried at 800°C for 10 minutes and immediately welded. Sealed 
capsules were packed with hematite powder into 3/4 inch 
furnace assemblies with tapered graphite heaters, which 
reduce the temperature gradient along the 10 mm capsule to 
15°C (Kushiro, 1976). The hematite acts as a trap for any 
water diffusing into the charge from the assembly. The 
assemblies were stored in a drying oven at 110°C prior to 
use. 

Temperatures were monitored with a Pt/Pt13Rh 
thermocouple without any correction for pressure, and are 
believed accurate to better than +10°C. Pressures were 
monitored continuously with a Heise gauge and were accurate 
to within +0.5 kbars. Pressure calibrations (by DTA) were 
performed using the melting curve of NaCl (Clark, 1959) and 
a pressure correction of -7% was applied to all runs. 

Run durations were one hour with the exception of a 
zero-time experiment and one 1/2 hour experiment (run 
numbers 12 and 10 (Table 2), respectively). The zero-time 
experiment provided confirmation of an initrablyttiat 
interface. Runs were quenched by switching off the power to 
the heater resulting in quench rates of approximately 


125°C/sec. Quenched runs (encased in coarsely recrystallized 
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specular hematite) were set in epoxy and sliced in half 
longitudinally. The platinum capsules preserved their 
cylindrical form with only minor necking during the run. 
Thin sections of the charges were examined optically 
revealing colorless, transparent glasses with no crystals 
(even for one run conducted below the liquidus of jadeite at 
1200°C and 15 kbars). The melt couple interfaces were 
optically invisible. 

The charges were analyzed after each experiment for Na, 
Al and Si by the energy dispersive technique described above 
and for fluorine by wavelength dispersive analysis. Fluorine 
analyses were standardized to a sample of the original 
fluorine-bearing jadeite glass. This glass had been analyzed 
against NBS opal glass SRM #91 by neutron activation 
analysis (chapter 1). The use of a fluorine standard with. 
the same matrix composition as the analyzed samples 
minimizes errors associated with poor ZAF correction factors 
for fluorine. 

Wavelength dispersive analysis for fluorine required a 
sample current of 40 nA to achieve reasonable count rates 
(approximately 40 cps/wt% F). Therefore, in order to avoid 
Significant volatilization during fluorine analysis two 
techniques were employed. In the first case, a point beam 
waS moved continuously, perpendicular to the profile, 
covering a distance of 50 micrometers in 100 seconds. In the 
second case, the beam was rastered over a 10 x 10 micrometer 


area for 100 seconds. The fluorine totals normalized to the 
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Standard were identical for both techniques and chart 
recordings showed that no significant volatilization of 
fluorine occurred during the Stationary analyses. Fluorine 
(WDA) and Na, Al and Si (EDA) spectra were combined and 


reduced using EDATA2 (Smith and Gold, 1979). 


C. Results 

This study involves chemical diffusion. Therefore, 
before the data are discussed and comparisons are made, we 
must distinguish between three categories of diffusion, 
namely, tracer, self and chemical diffusion. Self-diffusion 
is the diffusion of a single chemical species in the absence 
of a chemical gradient. This type of diffusion is usually 
investigated by labelling some of the diffusing atoms with 
an isotopic tracer (e.g. Shimizu and Kushiro, 1984). Tracer 
diffusion is the diffusion of an individual species in the 
presence of a small but finite chemical gradient. The tracer 
is usually a radioisotope (e.g. Watson, 1979) and the tracer 
concentration is so small that no significant concentration 
gradients result for any of the other melt components. 
Chemical diffusion is the diffusion of two or more species 
in response to a chemical activity gradient (e.g. Kushiro, 
1983). Tracer and self-diffusion are practically equivalent 
for natural melts and usually represent minimum 
diffusivities for the components being studied because there 


is no large chemical activity gradient to serve as a driving 
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force for diffusion. The similarity of tracer and 
self-diffusion is evidenced by the fact that they are both 
described mathematically by the single component form of 
Fick's first and second laws (Crank, 1975; Hofmann, 1980). 
In a general sense, chemical diffusion is a multicomponent 
phenomenon, but it may be approximated to a binary 
interdiffusion process when the chemical activity gradient 
for all other species (and the resulting diffusion of these 
Species) is insignificant. For such cases, an effective 
binary diffusion coefficient (EBDC) may be obtained (Cooper, 
1968) and we may talk in terms of a binary interdiffusion 
coefficient. 

In the present study, quantitative analyses for Na, Al 
and Si revealed no concentration gradients for these 
elements, indicating that the diffusion process could be 
approximated by binary interdiffusion of fluorine and 
oxygen. ‘Therefore, diffusion profiles of fluorine (Fig. 9} 


were fitted to the following form of equation (Crank, 1975): 


] 


X =ferf ([2C - (C,+C) 
au ao 


where D is the interdiffusion coefficient (cm’*/sec) , t is 
time (seconds) and C, C; and C2 are the concentration at 
distance x, the maximum and the minimum concentrations of 
fluorine, respectively. Erf~' is the inverse of the error 
function. The interface (x=0) was optically invisible and, 


therefore, the half maximum of the Giffusion profile was 
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Figure 9. Diffusion profile of fluorine in jadeite melt. 
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chosen to represent the interface. This choice was confirmed 
by the fit of the error function (Fig. 10). Identical 
diffusion profiles at several locations across each charge 
ruled out the possibility of significant deformation of the 
interface during individual runs. Plotting 

Brome C- (C7 7*C2) 1) /(C."Cl)). ws. ax Rig a 10 ivields linear 
plots whose slopes atan/o(be. A linear dependence of 
erf~-'((C-(C,+C2))/(C,-C2)) on x indicates that the diffusion 
process is independent of concentration. The values of D 
obtained from each run are presented in Table 4. 

The data of Figure 11 show increasing diffusivity of 
fluorine with increasing temperature. If we assume a linear 
dependence of log D on reciprocal temperature, we may fit 
the data of Table 4 to an Arrhenius equation for each 
pressure. The temperature dependence of diffusion is 


represented by the following form of equation: 


Tog, 9D = 10g, Dy - Ea/2.303RT 


where D is the diffusion coefficient at temperature T (K), 
D,o is the Arrhenius frequency factor, Ris the gas constant 
and Ea is termed the Arrhenian activation energy of 
diffusion. The results of least squares fits to the data in 
Table 4 are presented in Table 5. F-O interdiffusion 
activation energies range from 36 to 39 kcal/mole. 

The pressure dependence of diffusion may be linearly 


approximated by an Arrhenius equation of the 
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Figure 10. Reduced diffusion profile using equation (1) and 
data from fluorine-rich (triangles) and fluorine poor 
(inverted triangles),limbs of the diffusion profile. The 
slope corresponds to/2{Dt (Data from run no. 3). 
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Figure 11. Diffusion data for 10, 12.5 and 15 kbars and 
computed Arrhenius parameters. 
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Table 4. Experimental conditions and results 


un Temperature Pressure Duration log D. _ # of points 


oe 

(°C) (kbars) (sec) Comey ec) 
] 1400 15 3600 -6.18+.04 16 
2 1200 15 3600 -6 .87+.03 25 
3 1300 15 3600 -6 .39+.01 33 
4 1400 © 10 3600 ~6.15+.03 20 
5 1200 10 3600 -6 .87+.02 14 
6 1300 10 3600 6.4% .02 18 
7 1300 12.5 3600 -6 .60+ .04 17 
8 1400 1245 3600 -6.18+.01 9 
9 1200 25 3600 -§ .87+.05 9 
10 1300 10 1800 -6.41+.06 1g 
fl 1400 10 3600 -6.27+4.02 il 
a2 1300 10 0 - - 
13 1300 15 3600 -6 .62+.03 10 


uncertainties in log D are quoted at 1 standard deviation. 
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Table5 . Arrhenius parameters 


Pressure Ea log Dy # of points 
(kbars) (kcal/mole) 


10 36+2 -1 .464+.3 ct 
1205 39+2 -1.18+.2 3 
15 39+4 -1.12+.6 4 
Temperature Va log Do # of points 
(2G) (cm? /mole) 
1200 On6 -6 .87 S) 
1300 -1.4+.9 -6.36+.09 5 
1400 0.941 -6.274.09 4 


Eee 
uncertainties quoted at 1 standard deviation 
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form: 


10g, 9D = 10g, Dp - VaP/2.303RT 


where P is the pressure (dyne/cm?) and Va is the 
Arrhenius activation volume (cm*/mole). Calculated values of 
Va obtained from the least squares fits to the data for 
1200, 1300 and 1400°C are presented in Table 5. The 
activation volumes are small, with large uncertainties, 
reflecting the extremely small pressure dependence of F-O 
interdiffusion. Therefore, a mean activation volume of -0.33 
cm?/mole was calculated (based on all twelve data points) 


assuming temperature-independence of the activation volume. 


D. Discussion 


Comparison with oxygen diffusion 

Shimizu and Kushiro (1984) have measured oxygen 
self-diffusivity in jadeite melt from 5 to 15 kbars and 1400 
to 1610°C. They have reported self-diffusivities ranging 
from 6.87 x 10°'® to 4.72 x 10°'° cm?/sec. These values are 
three orders of magnitude lower than the F-O 
interdiffusivities at the same pressure and temperature and 
this difference shows that the presence of anionic chemical 


activity gradients can result in a large increase in oxygen 


diffusivity. 
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Comparison with tracer diffusion 

Tracer diffusion of fluorine has been studied by 
Johnston et a]. (1974) in a eutectic composition in 
CaO-Al1203-Si02, that contains approximately 40wt% CaO, 20wty% 
Al203 and 40wt% SiO2, (CagoAle2oSigg). The results of this 
study are useful for a discussion of fluorine diffusivity 
because CayoAlzoSiugo represents a depolymerized melt which 
has been well-studied in the glass literature, and the 
diffusivities of Ca, Al, Si, and O have been measured. 
Figure 12 presents the results of tracer diffusion studies 
on Ca (Towers and Chipman, 1957), Al (Henderson et ai/., 
1961), Si (Towers and Chipman, 1957), and O (Koros and King, 
WoG2mmOishi ef a@]/., 1975))-in CagoAlsoSiso amelt.) lb as clear 
from Figure 12 that fluorine diffuses faster than any of the 
other species studied. Also, fluorine tracer diffusivity has 
the lowest activation energy of any of the elements studied. 
As stated above, chemical diffusivities are usually larger 
than tracer diffusivities; however, the magnitude of tracer 
diffusivity in CayoAlzoSiyo melt is remarkably large. The 
observation that tracer diffusion in CayoAlzoSiuo melt is 
faster than chemical diffusion in jadeite melt implies a 


strong composition dependence of fluorine diffusivity. 


Comparison with Si-Ge and Al-Ga interdiffusion 
Kushiro (1983) has investigated Si-Ge and Al-Ga 


interdiffusion in jadeite melt from 6 to 20 kbars at T2005G. 
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Figure 12. Tracer diffusivities of various ions in 
lime-aluminosilicate melt (see text for data sources; 
Ben. oxygen diffusivity data is from Koros and King, 
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A marked asymmetry is observed in Si-Ge and Al-Ga 
interdiffusion profiles indicating a strong dependence of 
Si-Ge and Al-Ga interdiffusivities on Si/(Sit+tGe) and 
Al/(Al+Ga), respectively (Kushiro, 1983). Comparison of 
these interdiffusion data with the data for F-O 
interdiffusion at 1400°C and 15 kbars shows, as expected, 
that F-O interdiffusion is much faster than either Si-Ge or 
Al-Ga interdiffusion. 

Kushiro (1983) did not investigate the temperature 
dependence of cationic interdiffusion, but he showed that 
the pressure dependence of Si-Ge interdiffusion was much 
larger than the pressure dependence of Al-Ga interdiffusion. 
Figure 13 is a plot of the Arrhenius activation volume, Va 
(cm?/mole), versus the ratio of the ionic volumes of the 
interdiffusing species. The ionic volumes were calculated 
using ionic radii from Whittaker and Muntus (1970) for O, F, 
Al, Si, Ga, and Ge. The Va data are least squares fits to 
he 400°C, data from Kushiro (1983) for Si{GerandeAl-Ga 
interdiffusion and the mean value of Va for F-O 
interdiffusion from this study. It iS apparent in Figure 13 
that the relative sizes of the interdiffusing species are a 
large factor in determining the pressure dependence of 
interdiffusion. Therefore, it is anticipated that the 
pressure dependence of interdiffusion Containsme tole m1: 


any, information on the pressure dependence of 


self-diffusion. 
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Figure 13. Pressure dependence of F-O, Al-Ga and Si-Ge 
interdiffusion. Activation volume (Va) versus the ratio 
of ionic volumes of the interdiffusing species (symbol 
size corresponds to 1 standard deviation uncertainty in 
the Va data for F-O interdiffusion). Inset. The relative 
pressure dependence of F-O, Al-Ga and Si-Ge 


interdiffusion (Al-Ga and Si-Ge data from Kushiro, 
1983). ! 7 
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Compensation relationships 

Figure 14 is a plot of the frequency factor versus the 
activation energy of diffusion. Figure 14 will be used for 
two purposes; firstly, for a discussion of the correlation 
between frequency factors and activation energies (i.e. 
compensation) and secondly, to illustrate the effects of 
fluorine and water on cationic diffusivities in silicate 
melts. For the discussion of compensation reference will be 
made to the lines labelled 1 through 5 and the data for F-O 
interdiffusion (open circles) and H20 chemical diffusion 
(solid circles). The inset and cationic diffusivity data in 
Figure 14 will be discussed in the next section. 

Winchell and Norman (1969) first showed a positive 
correlation between the frequency factors and activation 
energies of various cationic diffusivities in CayoAl2zoSiyo 
melt (Fig. 14; line 2). From this correlation they defined a 
compensation "law" for silicate melt diffusivities which 
Winchell (1969) extended to other synthetic melt 
compositions. Hofmann (1980) proposed an equivalent 
relationship for basalt and obsidian melts incorporating 
cationic and oxygen diffusivity data (line 4). Hart (1981) 
distinguished two separate compensation laws, one for basalt 
melts (line 1), and one for obsidian melts (line 5), and 
finally, Dunn (1982) proposed a compensation law for various 
Synthetic silicate melts based on oxygen diffusivities (line 
3). The most important physical implication of compensation 


ies viicate Liquids 1s that all melt diffusivities converge 
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Figure 14. Compensation plot for silicate melt diffusivities 
showing 1) data for H20 chemical diffusion (Karsten et 
al., 1982: solid circles)and F-O interdiffusion (this 
study; open circles) 2) compensation laws (solid lines) 
3) cationic diffusivities in dry, fluorine free melts 
(squares); in hydrous melts (triangles); and in 
fluorine-bearing melts (diamonds). 
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at a unique temperature (T critical or Tc), which may or may 
not be attained in nature. The implications of such behavior 
are discussed at length by Hart (1981). Also, compensation 
plots provide a test of the consistency of diffusion data 
for a given melt composition. Figure 14 is a compensation 
plot which includes the results of this study for the 
chemical diffusion of fluorine and the results of Karsten et 
al. (1982) and Delaney and Karsten (1981) for the chemical 
diffusion of water. It is clear in Figure 14 that the 
fluorine data of this study may be included equally well 
within any of the compensation law relationships except the 
Obsidian curve of Hart (1981). The same conclusion is 
evident for the H20 data. The agreement of the fluorine and 
water data with the compensation laws for depolymerized 
melts (lines 1 to 4) suggests that water and fluorine 


depolymerize obsidian and jadeite melts, respectively. 


Effect of water and fluorine on cationic diffusivities 
Figure 14 emphasizes a second and very important aspect 
of fluorine and water in silicate melts. Included in Figure 
14 are data on the effects of water on Ca, Na and Cs 
diffusivities (Watson, 1981) and fluorine on Ca and Fe 
diffusivities (Johnston et a/]., 1974). The influences of 
water and fluorine on cationic diffusivities are large. All 
cationic diffusivities investigated increase when either 


water or fluorine is added to the Silicate melt. In most 
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cases, the increased diffusivity yields a lower activation 
energy and frequency factor (log; oD.) and, in these cases, 
the Arrhenius relationship for cationic diffusivity 
describes a line which rotates about a critical temperature, 
Tc (Fig. 14, inset, case A). This behavior produces a trend 
on the compensation plot which is sub-parallel to the 
various compensation relationships. In Figure 14, several 
cationic diffusivities are affected in this manner (e.g. 
increases in Ca and Fe diffusivities in fluorine-bearing 
melt and increases in Cs and Na diffusivities in hydrous 
melts). The behavior of Ca in hydrous melt is somewhat 
anomalous because the increase in Ca diffusivity takes the 
form of a bulk translation of the Arrhenius line (Fig. 14, 
inset, case B) yielding a lower activation energy but a 
higher frequency factor. However, the trends of Ca and H20 
diffusivity with increasing water content are similar (Fig. 
14). Regardless of the mechanisms responsible, increased 
diffusivity of cations with addition of fluorine or water to 
Silicate melts has significant implications for the roles of 


fluorine and H,O in establishing chemical equilibrium during 


igneous processes. 


Geological applications 
Jadeite melt has a highly polymerized structure Similar 
to natural melts whose compositions are approximated by the 


system nepheline-kalsilite-silica (Seifert et al]., 1982). 
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Relatively dry, rhyolitic, trachytic and phonolitic melts 
often contain up to 1 wt% fluorine (Carmichael et al .7e@ 1974; 
Bailey; °1977;° Christiansen et al.,' 1983). The present 
diffusion data provide information on a very important 
aspect of the petrogenesis of relatively dry, felsic melts. 
It has been observed experimentally that dry or 
water-undersaturated melts of rhyolitic and feldspathic 
composition have extremely slow equilibration rates due to 
low diffusivities in the melts (Schairer, 1950; Piwinskii, 
1967; Whitney, 1975). Johannes (1978, 1980) has proposed 
that even water-saturated granitic melts have equilibration 
rates, below 700°C, which are low enough to yield metastable 
melt compositions in nature. Considering the dramatic effect 
of fluorine on diffusivities in silicate melts, the presence 
of fluorine in dry rhyolitic melts could be a crucial factor 
in determining the rate of establishment and the physical 
extent of chemical equilibrium during anatexis in the lower 
Ceust . 

Fluorine increases cationic diffusivities in 
depolymerized melts (such as Cayo AlSOS140%" Johnston ets a] 15 
1974). Therefore, perhaps fluorine is capable of increasing 
diffusivities in melts which are already depolymerized due 
to dissolved water. This potential additive effect of 
fluorine and water on the diffusivities of various cations 


in late-stage, water-saturated, granitic melts may yield 


exceptionally high melt diffusivities. 
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E. Summary and conclusions 

The chemical diffusion of fluorine in a jadeite melt 
involves binary interdiffusion of fluorine and oxygen. This 
interdiffusion is concentration independent, in contrast to 
Si-Ge and Al-Ga interdiffusion in jadeite melt. High 
temperature fluorine diffusion, both chemical and tracer, is 
equal to, or greater than both cationic and oxygen diffusion 
in fluorine-free melts. The chemical diffusivity of fluorine 
in jadeite melt is similar in magnitude to the chemical 
diffusion of water in obsidian melts: however, melt 
compositional dependence of fluorine diffusivity precludes a 
direct comparison of water and fluorine diffusivities at 
this time. The chemical diffusion of fluorine has an 
Arrhenius activation energy of 36-39 kcal/mole compared with 
19 kcal/mole for chemical diffusion of water in obsidian 
melts and for tracer diffusion of fluorine in CayoAl2ocSiuno 
melt . The results fit several of the compensation "laws" 
which have been proposed for quantification of cationic and 
anionic diffusivities in depolymerized silicate melts. A 
Significant effect of fluorine on melt diffusivities is that 


cationic diffusivities are enhanced by the addition of 


fluorine to silicate melts. This behavior is also observed 


when water is added to silicate melts. The effect of 


fluorine on melt diffusivities may be a important factor in 


the chemical equilibration of dry, igneous melts. 
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IV. One atmosphere diffusion 


A. Introduction 

In chapter 2, the comparison of fluorine diffusivity in 
the lime-aluminosilicate melt and jadeite melt indicated 
that the value of fluorine diffusivity may be strongly 
dependent on melt composition. The experiments in this 
chapter were conducted to 1) investigate the melt 
compositional dependence of F-O interdiffusivity, 2) further 
investigate the pressure dependence of F-O interdiffusion in 
jadeite melt and 3) compare the chemical diffusivities of 
water in obsidian and fluorine in albite auaice pressure. 

The volatilization of fluorine from silicate melts has 
been investigated by several workers uSing weight loss and 
bulk chemical analysis techniques (Kumar et a/., 1961; 
Barlow, 1965; Kogarko et a]., 1968; Al-Dulaimy, 1978). The 
present study uses the relatively high volatility of 
fluorine in melts in the system Na20-Al203;-Si02 to 


investigate fluorine diffusion. 


B. Experimental method 


The starting materials for this study were the glasses 
used for the viscometry work in chapter 1. Chips of glass 
were ground into spheres ( 2-4 mm in size) using the 
technique of Bond (1951). Spheres of glass were suspended 


(with white glue) from platinum loops iMitan 


electrically-heated vertical tube furnace, equipped with a 
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gas-tight alumina muffle tube. An "infinite" O2 gas 
reservoir was maintained during the experiments by flowing 
oxygen gas through the furnace at a linear flow rate of 0.10 
cm/sec. Temperatures were controlled using a Pt-Pt13Rh 
thermocouple which was suspended close to the spheres. Three 
spheres were run at each temperature for each composition. 
The experiments ranged in duration from 1200 to 19,200 sec. 
and melt spheres were quenched in air by removal from the 
furnace. Quenched melt (glass) spheres were removed from the 
ceramic hanger, ground in half, mounted in epoxy and 
polished for electron microprobe analysis. In rare cases 
where the glass beads deformed from spherical shape, due to 
improper loop size or mounting, the runs were discarded. 
Generally, spherical shape was extremely well-preserved. 
Some runs exhibited surface bubbles, presumably due to vapor 
formation near the melt-O, interface. These bubbles did not 
interfere significantly with the following analysis because 
the bubbles were restricted to a thin surface layer (see 
below). 

Analyses of the concentration profiles were performed 
using an ARL-SEMQ electron microprobe. Operating conditions 
for the scans of Na, Al and Si were a 15kV accelerating 
voltage and a 4 nA sample current with the beam rastered 


over a 10 by 10 micron area. Fluorine scans and step 


analyses required a sample current of 40 nA to achieve 


reasonable count rates (approx. 40cps/wt% F). The 


quantitative analyses for fluorine were 100 sec counts using 
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a rastered beam as for Na, Al and Si. 

Wavelength dispersive scans for Na, Al, Si and fluorine 
(Fig. 15) show two basic features. Firstly, the 
concentration profile of fluorine extends (over several 
hundred microns) from essentially zero concentration at the 
Sphere edge to a plateau corresponding to the original 
undepleted fluorine concentration of the melt. Secondly, the 
concentrations of Na, Al and Si are constant over almost the 
entire fluorine depletion zone except within 30-50 microns 
of the sphere edge. 

Quantitative analyses of fluorine content were obtained 
by standardizing the raw counts from step scans across each 
Sphere to the undepleted central zone of the sphere. The 
exchange of up to 6 wt.% fluorine for oxygen has no effect, 
within error, on the fluorine ZAF correction. Consequently, 
after background correction, raw counts from step traverses 
to the edge of each sphere could be normalized to the counts 
from the center of each sphere to obtain quantitative 
fluorine analyses. 

The 30-50 micron wide "surface layer" of Si depletion 
and NatAl enrichment strongly suggests the escape of 
fluorine from the melt surface is aS SiF, gas. This 
observation is corroborated by previous studies of fluorine 
volatility in dry systems at 1 atmosphere (Kumar et al., 


1961: Barlow, 1965; Al-Dulaimy, 1978). 
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Figure 15. Microprobe scans for Na, Al, Si and fluorine 
across sphere edge. 
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C. Data reduction 

The constant levels of Na, Al and Si inside the surface 
layer mean that bulk diffusion of fluorine to this layer may 
be treated as binary F-O interdiffusion. The mathematical 
treatment of diffusion for these experiments is for a 
one-dimensional semi-infinite medium (Crank, 1975; p.35-38). 
The chemical activity gradient which serves as the driving 
force for chemical diffusion is quantified by the "reduced 
concentration" term (C-C,)/(Co-C,) where C is the measured 
concentration along the profile, C, is the initial 
concentration in the melt, and Co is the concentration in 
the gas phase. This equation holds equally well for the case 
of absorption into a semi-infinite melt medium from a gas 
phase (e.g. the hydration experiments of Shaw, 1974) as it 
does for the present case of desorption or volatilization. 
In these volatilization experiments, Co = 0, because the gas 
phase is essentially pure O, gas. Therefore, the reduced 
concentration term simplifies to 1 - C/C,. In the case of 
concentration-independent diffusion, the equation relating 


the diffusion coefficient (D), distance (x), and time (t) is 


as follows: 


ener (6/£,) 


2tbt 


where erf-' is the inverse of the error function and x is 


the distance from the melt-vapor interface. In practice, 


plots of erf-‘(C/C;) vs. x yield straight lines whose slopes 
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equal 1/(2JDt). A typical plot is shown in Figure 16. 


D. Results 

The experimental results are summarized in Table 6 and 
plotted versus reciprocal temperature in Figure 17. The 
temperature dependence of F-O interdiffusion is fitted to 


Arrhenius equations of the form: 


Tog, 5D = 1091 QD, - Ea/2.303RT 


where D is the diffusivity at temperature T, Dy is the 
pre-exponential or "frequency factor", R is the gas constant 
and Ea is termed the activation energy of F-O 
interdiffusion. The Arrhenius parameters, log; oD. and Ea, 
for each melt composition are listed in Table 7. 

The interdiffusion of fluorine and oxygen in the system 
Na20-Al,03;-Si0O, is strongly dependent on melt composition. 
Ae v400°C,), F-O interdiffusion increases ‘over two orders ‘of 
magnitude in melts along the join NaAlO2-Si02z from albite 
(logioD = -8.46) to nepheline (logioD = -6.26). F-O 
interdiffusivity also varies strongly with alkali/aluminum 
ratio. At 1200°C and 75 mole % SiOz, the diffusivity of 
fluorine in albite is slower than in the peralkaline and 
peraluminous melts. 

The activation energy of F-O interdiffusion increases 
along the join NaAlO,-SiO2 from albite (29 kcal/mole) 


through jadeite (34.4 kcal/mole) to nepheline (42.4 
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Figure 16. Results of the time series experiments on jadeite 
melt at 1200°C. 
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Figure 17. Results of the 1 atm diffusion experiments 
(asterisk denotes experiments conducted below the 
liquidus of the fluorine-free melts; square symbol 
denotes data for the peralkaline melt). 
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Table 6 : Experimental conditions and results 


Composition Temperature Duration -logaeb # of expts. 
(Celsius) (min) (cm? /sec) 

Jadeite 1200 10-320 Ee 7 
1300 3600 7224 3 
1400 1860 ren 3 
Albite 1200 6900 8.98 3 
1300 7680 8.84 3 
1400 7200 8.46 3 
Nepheline 1200 1920 He02 3 
1300 1380 6.74 S 
; 1400 1140 Gaz 6 3 
Peraluminous 120 11280 O50 3 
1300 7560 7.98 3 
1400 7680 7.74 3 
Peralkaline 1200 4260 7.40 3 


x Results of a time series of seven experiments (see figure 

he we-uncertainty in -1lo0g;9D is estimated to, bevt0.\slog 

Upiesmbasea onsthei stamdardydevaationsof. the seven time 
series experiments in jadeite at 1200°C. 


Table 7 : Arrhenius parameters for- li atmosphesen EO 
interdiffusion 


Tae i D Ea 
Composition Feces NEES (kcal/mole) 
Albite meee sn 
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kcal/mole). If the trend of decreasing diffusivity with SiO, 
content extends beyond albite along the join NaAl0O2-Si0>2 
then the F-O interdiffusivity measured for albite may be 
taken as an upper limit on F-O interdiffusivity in obsidian 


melts. 
E. Discussion 


Comparison with viscosity data 

The suggestion that similar mechanisms and species may 
control the viscous flow and diffusivities in silicate melts 
has led several investigators to examine the validity of the 
Stokes-Einstein equation which inversely relates viscosity 
and diffusivity (Watson, 1979; Hofmann, 1980; Jambon, 1982; 


Burnham, 1983; Shimizu and Kushiro, 1984). It is: 


Dn = KbT/6yxr 


where Kb is the Boltzmann constant and r is the radius of 
the diffusing species. This form of equation predicts a 
linear inverse correlation between log D and log n. It has 


been generally concluded that the Stokes-Einstein 
relationship is qualitatively invalid in relating viscosity 


and network-modifying cationic diffusivities in Silicate 


melts (Winchell and Norman, 1969; Magaritz and Hofmann, 


1978: Watson, 1979; Hofmann, 1980; Jambon, 1982; Angell et 


al., 1983). However, for the case of oxygen diffusion, 
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and Kushiro. (1984) have stated that an inverse correlation 
with viscosity is well approximated by the Eyring form of 
the Stokes-Einstein equation. The success of any equation 
which inversely relates viscosity and oxygen diffusivity 
Strongly implies that the same structural unit and mechanism 
are involved in both transport processes. Dunn (1983) has 
recently shown that 0?” is probably the dominant species 
involved in oxygen diffusion in basaltic melts and Shimizu 
and Kushiro (1984) argue that the viscous flow of jadeite 
and diopside melts is controlled by the diffusion of 
individual 0?” ions. 

The conclusions of this study regarding the application 
of the Stokes-Einstein equation to F-O interdiffusivity are 
negative. Correlation of log n versus log D 1S very poor. 
Qualitatively, we would expect this result because the 
viscosity of these melts is a strong function of fluorine 
content but the F-O interdiffusivity is not. There is even 
an increase in F-O interdiffusivity with increasing 
viscosity from the albite to the peraluminous melt. Clearly 
the presence of fluorine in these melts removes the 


applicability of the Stokes-Einstein equation to the 


diffusion of oxygen. 


Pressure dependence of F-O interdiffusion in jadeite melt 
The results of this study for jadeite melt are compared 


with the high pressure data of chapter 2--in Pigure W6sen-0 


interdiffusivity clearly increases with pressure from 0.001 
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Figure 18. Pressure dependence of fluorine diffusion in 
jadeite melt.Inset: Pressure dependence of activation 
energy for F-O interdiffusion. 
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to 10 kbars, as does oxygen self-diffusion (Shimizu and 
Kushiro, 1984). The lines are interpolations based on the 
assumption of a smoothly decreasing pressure dependence with 
increasing pressure. In fact, Dunn (1983) has shown that the 
pressure dependence of oxygen chemical diffusivity in 
basaltic melts is discontinuous in the pressure range of 1 
to 10 kbars whereas, Shimizu and Kushiro (1984) show a 
smoothly increasing self-diffusivity of oxygen in jadeite 
-melt in the pressure range of 5 to 20 kbars. This contrast 
may result, in part, from comparison of chemical and self 
diffusivities and thus the pressure dependence of F-O 
interdiffusion in jadeite melt is difficult to discuss a 


structural terms. 


Comparison with water 

The comparison of the chemical diffusion of fluorine 
and water is-interesting for twosreasons..Firstly,,the 
geological significance of fluorine in natural melts was 
discussed in chapter 2. It is clear that relatively dry, 
fluorine-rich geological melts do exist and their kinetic 
behavior requires investigation and comparison with that of 
wet melts. Secondly, the structural comparison of fluorine 
and water in chapter 1 led to the conclusion that the 
depolymerizing influences of fluorine and water have very 
similar effects on the viscosity of silicate melts. 


Chemical diffusion of water in natural silicate melts 


has been investigated by several workers (Shaw, 1974; 
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Friedman and Long, 1976; Arzi, 1978; Jambon et al]., 1978; 
Delaney and Karsten, 1981; Karsten et a]., 1982): In Figure 
19, the data for the chemical diffusion of fluorine in 
albite melt are compared with the available data for the 
chemical diffusion of water in obsidian melts. Although the 
comparison of data in Figure 19 involves melts of albite and 
obsidian composition both represent relatively polymerized 
melts with alkali/aluminum (molar) ratios at or near 1:1. 
The comparison of Figure 19 uses 1 atmosphere data for 
Fluorine and 0.1 - 2.0 kbar data for water diffusion. With 
these provisions in mind, this comparison is discussed 
below. 

The chemical diffusivity of fluorine in albite melt is 
approximately two orders of magnitude less than that of 
water in obsidian melt. Considering the composition 
dependence of fluorine diffusivity observed from nepheline 
melt tio albite melt (Fig. 17), the difference™between 
fluorine and water diffusivities in obsidian is probably 
even larger. The temperature dependence of fluorine 
diffusivity yields an activation energy of 29 kcal/mole in 
albite melt compared with an activation energy of 19 
kcal/mole for water diffusion in obsidian (Delaney and 


Karsten, 1961). Such differences in diffusivity and 


activation energy may result from different mechanisms of 


a@iffusion for water and fluorine in these polymerized melts. 


Qualitatively, these differences imply that tier so ecne 


principal diffusing species, as is suggested by the binary 
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Figure 19. Chemical diffusivities of fluorine in albite melt 
and H2O0 in rhyolitic melts. (see text for data sources) 
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Mature of interdiffusion of fluorine and oxygen, then water 
1s probably not transported principally as OH-, the 
monovalent anion of similar size, Conversely, if the 
diffusion of water is as OH-, then the species involved in 
fluorine diffusion might be a larger complex .such as AlF,*- 
(Manning et a/]., 1980). 

Perhaps the most significant difference between the 
diffusive behavior of fluorine and water is the 
concentration dependence of the chemical diffusion of water 
in obsidian melt contrasted with the concentration 
independence of fluorine diffusion in all melts investigated 
in this study. The concentration dependence of water 
diffusivity in obsidian melt has been most recently 
investigated by Delaney and Karsten (1981) and Karsten et 
al. (1982). These workers state that the activation energy 
of chemical diffusion of water remains constant while the 
frequency factor (log D.) increases with water 
concentration. This increase has been interpreted (Karsten 
et al., 1982) to result from an increase in the average jump 
distance associated with the diffusive mechanism. Karsten eft 
al. (1982) further propose that the concentration dependence 
of water diffusion results from the occurrence of at least 
two distinct solution sites for water in obsidian melt whose 
relative occupancies change with water concentration. 

In addition, the concentration dependence of water 


diffusivity, combined with the effect of water on melt 


viscosity, permits application of the Stokes-Einstein 
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equation to melts of constant composition and varying water 
contents (Burnham, 1983). The concentration independence of 
fluorine diffusion contrasts strongly with the behavior of - 
water and suggests that a single structural site is entirely 


.adequate to explain the observed diffusive behavior. 


Diffusion mechanisms 

From a combination of data of viscosities and 
diffusivities in F- and H,0-bearing melts (chapters 1,2 and 
above) it is apparent that the effects of fluorine and water 
on melt viscosity are similar, whereas their diffusivities 
are considerably different. The viscometry data implies 
similar structural roles for fluorine and water within 
polymerized silicate melts, probably the replacement of 
Si-O-(Si,Al) bridges with (Si,Al)-F and (S1,Al)-OH, 
respectively. 

The diffusivity data, in contrast, highlight several 
differences in the transport of fluorine and water in these 
melts. Fluorine exchanges with oxygen via a mechanism which 
is independent of concentration and strongly dependent on 
melt composition. The exchange is effectively binary, 
probably without the involvement of any cations. 

Water, in contrast, diffuses in polymerized melts at a 
rate which varies with water concentration. Karsten ef al. 
(1982) have reported preliminary evidence of K and possibly 


Na concentration gradients produced during chemical 


diffusion of water. These alkali concentration gradients 


15 


suggest that interaction of water and alkalies may be an 
essential characteristic of chemical diffusion of water in 
these melts. The albite-H,0 water solubility model of 
Burnham (1975, 1979) does predict alkali transport during 
the solution of water in albite melt. 

Figure 20 presents tracer diffusion data for Li, Na, K, 
Rb and Cs (Jambon, 1982) and chemical diffusion data for 
water (Fig. 6 in Hofmann, 1980) in obsidian melts, as well 
as chemical diffusion data for fluorine in albite melt (this 
study). The diffusivity of alkalies in obsidian melt 
increases smoothly with decreasing ionic radius from Cs to 
Na. However, Li diffuses at essentially the same rate as Na. 
Jambon (1982) has interpreted the equivalence of Na and Li 
Srefusivity to result from the control, of Li*difttusivityeby 
the structure and mobility of the matrix or neighbouring 
cations with which Li will interdiffuse. Thus, in obsidian 
melt, the diffusion of Li may be controlled by the mobility 
of Na (Jambon, 1982). In Figure 20, it is evident that the 
chemical diffusivity of water, unlike that of fluorine, is 
very similar to the tracer diffusivity of Na (and Li). 
Perhaps the diffusion of hydrogen-bearing species in 
obsidian is controlled by the mobility of Na via an alkali 


exchange reaction similar to that described by Burnham 


(1975) for water solubility. 
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Figure 20. Alkali tracer diffusivities in obsidian melt 


compared with chemical diffusivities of fluorine and 
water. 
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V. Water solubility 


A. Introduction 

The solubility of water in granitic liquids of various 
compositions is a crucial property influencing the 
petrogenesis of intrusive and extrusive silicic rocks (e.g. 
Wyllie, 1979; Burnham, 1979). The dependence of solubility 
on melt composition is especially important during 
crystallization of granitic intrusions. Saturation of 
residual melts results in boiling phenomena closely 
associated with Cu-Mo porphyry systems and various 
weaneuet es kThe explosive eruption of magmas may be 
triggered by energy released during the expansion of a vapor 
phase evolved from melts (e.g. Sparks et @/]., 1977; Harris, 
1981b). Numerical models for the evolution of vapor (e.g. 
Wilson et a/., 1980) from granitic liquids at low pressures 
require detailed knowledge concerning the effects of 
composition and pressure on the solubility of water in such 
melts. 

The solubility of water in granitic melts has been the 
topic of several experimental studies; however, there is no 
consensus on the results or their interpretation (e.g. 
Oxtoby and Hamilton, 1978b; Day and Fenn, 1982). These 
previous solubility determinations may be divided broadly 
into two categories: those using chemographic or phase 
equilibrium techniques (Burnham and Jahns, 1962; Fenn, 1973; 


Whitney, 1975; Vogt, ef al., 1981) and those using a weight 
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loss method (e.g. Goranson, 1931, 1938: Bowen and Tuttle, 
1950; 2voderseteal gped957+ -Orlovas +1963; Oxtoby and Hamilton, 
1978a). Solubility data obtained by such techniques were 
reviewed by Day and Fenn (1982) who discussed the 
assumptions and limitations associated with each method. A 
new capacitance manometric technique for micro-determination 
of H20 evolved from vacuum fusion of glass (Harris, 1981a) 
affords a method for measuring H:;0 released from 
experimentally quenched vapor-saturated melts prepared as 
glass wafers free of fluid-filled vesicles. By using this 
technique the solubility of water in haplogranitic melts in 
the system K20-Na,0-Al203;-Si0z has been determined. 

The other unique aspects of this work are : (1) the 
study of haplogranitic melts to bridge the gap between melts 
of feldspar composition and natural granitic liquids; (2) 
the precise measurement of water solubility at up to four 
pressures in the range 1-2 kbars and determination of the 
pressure dependence for two compositions; and (3) 
investigation of the effects of increasing peralkalinity on 
water solubility at 1 kbar. The effects of composition and 
pressure on the solubility of H20 are essential for 


understanding the petrogenesis of silicic plutonic and 


volcanic rocks. 
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B. Experimental and analytical procedures 

The glasses used in the hydrothermal runs were 
synthesized from reagent grade Na7zCO7)h7 CO.) Ale On eands 
SiOz by using standard techniques described by Schairer and 
Bowen (1956). The compositions of the anhydrous glasses used 
as starting materials (Fig. 21) were determined by electron 
microprobe. 

Synthetic vapor-saturated melts were produced 
hydrothermally. Twenty-five to fifty milligrams of powdered 
glass were added to 10 to 20 milligrams of triple-distilled 
water in 3 mm Pt capsules and welded. The Apatite were 
checked for possible leakage by testing for weight loss 
after drying in an oven at 110°C for 15 minutes. Sealed 
capsules were loaded into cold-seal pressure vessels and 
raised to the desired run pressure by uSing an air-driven 
hydraulic pump. Pressurized capsules were heated to the run 
temperature (800°C) by using Kanthal-wound resistance 
furnaces mounted concentrically about the vessels along a 
vertical axis. Temperatures were monitored with 
Chromel-Alumel thermocouples and the actual temperatures are 
believed accurate to +5°C. Pressures were measured with 
Bourdon-tube gauges previously calibrated with a Heise 
gauge. All run durations were 10° seconds. The charges were 
quenched isobarically by using compressed air. Pressure was 
maintained during the quench by opening the vessel to the 


pressure line and by pumping for the duration of the quench. 


Owing to the response time of the pump there were slight 
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Figure 21. Melt compositions in the system 
K20-Na20-Al203;-SiO2 (1 kbar phase equilibria data from 
Tuttle and Bowen (1958) and Hamilton and MacKenzie 
(1965)). 
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fluctuations in pressure during quenching. These 
fluctuations were monitored and are included in the pressure 
uncertainties listed in Table 8. Except for water, the 
compositions of the hydrated glasses were determined by 
electron microprobe (Table 8). All of the glass 
compositions, when recalculated H,0-free and normalized to 
100%, are within analytical error of their respective 
anhydrous starting compositions. 

The quenched run products consisted of vesicular glass 
free of microlites. The vesicles contained liquid water and 
an enclosed vapor bubble and indicate that the charges were 
fFluid-saturated during the experiments. It is extremely 
difficult to accurately estimate the volume fraction of 
vesicles in the run products. This difficulty was avoided by 
preparation of vesicle-free glass wafers. 

Polished glass wafers free of occluded fluids were 
prepared manually from the experimental charges. The 
resulting wafers were typically 40-100 micrometers thick, 
depending upon the sizes and locations of vesicles within 
the glass wafer. Absence of occluded fluid was established 
by transmitted light microscopy at 200X magnification. The 
glass wafers were liberated from the Canada balsam mounting 
media by using acetone. The wafers were washed 
ultrasonically in acetone and checked again optically for 


purity at 200X magnification using transmitted and reflected 


light microscopy. 
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The amount of H20 dissolved in the glass was determined 
by the vacuum fusion micromanometric method described in 
Harris (1981a). The determinations of H20 were made on 
duplicate glass wafers of different mass for each run and 
the amounts reported are the duplicate means. The samples of 
glass ranged in mass from 0.020 to 0.600 mg. Sample masses 
were determined with a Cahn-G electronic microbalance and 
are considered precise to +0.003 mg. The analytical 
precision for water can be estimated by considering the 
variance of the relative deviations of the individual 
measurements about their duplicate mean. This method 
provides an estimate of +2.8% of the concentration for the 
standard deviation of the duplicate means about their 
ensemble mean, if multiple pairs of determinations were made 
for each experimental charge. The overall reproducibility of 
these solubility measurements includes uncertainties due to 
variation of run pressures, possible variable diminution of 
dissolved water during the quench, and analytical errors in 
the determination of H20 (e.g. incomplete extraction of H20 
during vacuum fusion of samples, or loss of H20 due to 
adsorption by the walls of the vacuum system). The overall 
reproducibility was estimated by repetition of experiments. 
Two solubility determinations have been repeated (Table 8; 
EA1.2K at 970 bars, runs 80 and 38; EA1.2K at 1620 bars, 
runs 85 and 44). Both solubility determinations were within 
0.06 wt% HO of their ensemble means, or about 1.4 % of the 


concentration. The precision attained in these measurements 
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exceeds that reported by Harris (1981a) Owing to the greater 
precision of sample masses that were made possible in this 


study by using an electronic microbalance. 


C. Results 
The glass compositions and solubility determinations at 
various pressures are given in Table 8. The effects of 


pressure and composition are described next. 


Pressure dependence 

The effects of pressure on the solubility of water were 
determined for two melts with identical Na/(NatK) ratios 
(0.57) and similar Si/Al ratios but with differing (Nat+K)/Al 
ratios. Both liquid compositions were peralkaline. The 
experimentally determined solubilities, expressed as weight 
percent H2,0 in the hydrous liquid, are shown in Table 8. The 
solubilities were recast as mole fractions by using the 
method of Burnham (1979) and fit (Table 9) to linear 
equations as functions of the square root of the water 
fugacity (Burnham et a]., 1969). The root mean square 
deviations from linearity are less than 1.4% (Figure 22). 
The equimolar solubility of water is slightly larger for the 
more peralkaline liquid at both 970 and 1620 bars P(H20). 


Hence, the derivatives with respect to Square root of 


fugacity differ by 5%. 
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Figure 22. Effect of pressure on water solubility. 
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TABLE 9 


EFFECT OF PRESSURE ON THE SOLUBILITY OF H,O IN Two Liguips 


ooo 


Pressure fx.o fu.o'- xm eae 
Run No. bars bars bars'- meas. calc. Residual 


EA!.2K 970-1630 bars 


80 970 825 28572 .363 .356 — .007 
38 970 825 28.72 .356 356 .000 
48 1160 965 31.06 .386 385 — .001 
43 1450 1175 34.28 426 425 = 001 
&5 1620 1297; 36.01 443 .446 .003 
net 1630 1305 36,12 437 447 O10 


XP = 0.003 + 0.0123 (fy.o. bars)? * 
r.m.s. error of fit = 0.005 


EAK 
36 970 825 DSes2 Sie rua OU 
83 1620 1297 36.01 .465 .465 .000 


X™ = 0.0006 + 0.0129 (fu.0. bars) * 
r.m.s. error of fit = 0.0006 


Notes.—Solubility data and the origin were fit to linear equations by the method of least squares. 
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Effects of Na-K exchange 

The effects of exchanging Na and K in haplogranitic 
melts may reveal whether there is a notable effect of alkali 
cation size on the solubility of water. Two experiments 
(runs 37 and 61) conducted at 970 bars show that the 
increase of the Na/K ratio from 1.34 at the eutectic to 2.23 
caused no measurable change in solubility (maximum 0.14 
wt%). Further research should determine whether this is also 
true for other liquid compositions in the vicinity of the 


ternary eutectic. 


Dependence of solubility on the alkali/aluminum ratio 

The equimolar solubilities of water in four liquids 
that differ in alkali/aluminum ratio were determined at 970 
bars and are summarized in Table 10. For the range of 
Silica-saturated compositions investigated, the equimolar 
solubility at 970 bars water pressure increases with 
increasing alkali to aluminum ratio and the root mean square 
deviation from linearity is about 1.4%. However, the 
solubility of water in the peraluminous liquid (Table 8, Run 
58) does not lie on the extension of this line to the 
peraluminous field. This suggests the existence of minima in 
the solubility of H,0 at (Nat+K)/Al=1 for melts along joins 


of constant molar $iOz in anhydrous melts (Rig 22237. 
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Figure 23. Effect of the alkali/aluminum ratio on water 
solubility. 
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TABLE 10 


DEPENDENCE OF THE EQUIMOLAR SOLUBILITY OF WATER IN MELTS ON THE ATOMIC RaTiO OF TOTAL ALKALIS 
TO ALUMINUM AT 970 Bars H;O PRESSURE 


SiO; (Na+ K)/Al 
Run Series wt % molar > Cus 
61 NOR 77.92 0.97 308 
ayy EOR TES 0.97 am) 
80 EA1.2K 77.49 1.19 : 363 
38 EAI.2K (71358 p20 356 
36 EAK P58 1.46 .372 


X™ = 0.200 + 0.124 (Na + K)/Al 
r.m.s. error of fit = 0.0110 
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Solubility in a phonolite liquid 

The 970 bar solubility of HO in the phonolite liquid 
exceeds that of all the granitic melts investigated. The 
phonolite composition, which has alkali/aluminum equal to 1 
and very much less silica than those runs listed in Table 
10, does not satisfy the linear relation obtained in Table 


tO 


Effect of fluorine on solubility 

Three fluorine-bearing melts of EOR composition, whose 
synthesis is Rencuited in chapter 6, were analyzed for water 
in an effort to determine the effect of added fluorine on 
the solubility. of water in this melt. The results of these 
analyses are presented in Figure 24. It is clear that 
fluorine added to haplogranitic melts decreases the 
solubility of water in these melts. These results contrast 
with the chemographic data of Koster van Groos and Wyllie 
(1968) and Wyllie (1979) which indicate enhanced water 
solubility in fluorine-bearing feldspar melts. The observed 
decrease in water solubility from this study implies that 
fluorine is incorporated within the silicate melt structure 
in a way which reduces the number of sites available for 
water dissolution. 

If water is dissolved predominantly as hydroxyl ions 
coordinated by tetrahedrally-coordinated cations (Burnham, 


1979, 1981) then perhaps fluorine is simply substituting for 
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Figure 24. Effect of fluorine on water solubility. 
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OH” in the depolymerized melt structure. 


D. Discussion 


Comparison with earlier determinations 

The determinations of water solubility in haplogranitic 
liquids from this study are shown (Fig. 25) in relation to 
other determinations at low pressures. Although the results 
are generally similar to those of earlier workers, the 
effects of pressure and composition on the water | 
solubilities are evident in the results of this study and 
may help to explain some of the differences shown in Figure 
Zan 

All previous measurements of water solubility in 
granitic melts (except those of Khitarov et al., 1959, 1967) 
have employed one or both of the chemographic and 
weight-loss techniques. The chemographic technique involves 
the construction of isobaric T-X or rock-water pseudobinary 
phase diagrams (see Burnham and Jahns (1962) for 
discussion). In these diagrams the intersection of the 
water-saturated liquidus with the divariant L-L+V or L-L+C 
curve is determined. The intersection gives the 
concentration of water equal to the solubility limit. 
However, according to Day and Fenn (1982), water contents 
may be overestimated by as much as 1.0 wt% by using the 


L-L+V or 'dimple' technique. At pressures below two kbar the 
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Figure 25. Solubility determinations in natural and 
synthetic granites. (BDM, Beinn an Dubhaich granite, 
Oxtoby and Hamilton, 1978b; SMG, Stone Mountain granite, 
Goranson, 1931; EDG, El'Dzurtinskii granite, Khitarov ef 
al., 1959; HP, Harding pegmatite, Burnham and Jahns, 
beo2). Data trom this study are for T= 600-6, 


o This study 
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L-L+C technique requires the careful addition of small 
amounts of water to experimental charges and the execution 
of water-undersaturated experiments on melts for which 
equilibrium is difficult to achieve. 

In the weight-loss technique, water concentrations are 
obtained by comparison of the weight of H2O added to the 
charge before the experiment with the weight of the post-run 
hydrous glass after the charge has been pierced and heated 
to 110°C (Hamilton et a]., 1964). Because the charges are 
generally oversaturated with H:0, the run products, 
especially below 2 kbar, are vesicular and the amount of H,0 
in the vesicles must be subtracted from the total H20 
determined by weight loss (Burnham and Jahns, 1962). 
Alternatively, if the vesicles can be ruptured by stepwise 
heating without loss of any H20 from the glass (e.g., Oxtoby 
and Hamilton, 1978a), then it may be possible to determine 
the solubility from weight loss measurements during heating 
of vesicular glass. Uncertainties for weight loss techniques 
are predictably rather large. 

In contrast, the experimental, sample preparation, and 
microanalytical techniques used in the present work permit 
the solubility to be determined with a variation of about 
+0.10 wt% (2 standard deviations) in replicate experiments 
at pressures from 1.0-2.0 kbars. The principal advantage 
over the previously mentioned techniques are the high 
precision and the fact that the measured water is wholly 


derived from the glass. Other significant advantages are 
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that the technique does not require multiple runs, precise 
determination of a phase boundary, or precise measurement of 
the amount of water added to an experimental charge. 

From the comparison of the present data and previous 
work for granitic melts depicted in Figure 25, two major 
conclusions can be drawn. The first conclusion is that the 
results of this study at 800°C and those of Khitarov et aij. . 
(1959) for the El'Dzurtinski granite at 900°C are 
approximately similar. On the other hand the data of Oxtoby 
and Hamilton (1978b) on the Beinn an Dubhaich granite, 
although similar in magnitude, were obtained at 1100°C. The 
second conclusion is that the solubilities from this study 
are lower than those of Burnham and Jahns (1962) on the 
Harding pegmatite at temperatures between 660-810°C and 
those of Goranson (1931) on the Stone Mountain at 900°C. The 
Stone Mountain, El'Dzurtinskii and Beinn an Dubhaich 
granites are approximately similar in alkali/aluminum ratios 
but they vary considerably in their K/Na ratios. The Harding 
pegmatite is distinct in having high contents of Rb, Li and 
fluorine (see Burnham, 1979). Oxtoby and Hamilton (1978b; 
see also Voigt et aJ., 1981) showed that the equimolar 
solubility of water is significantly greater in NaAl1Si30. 


melt than in KA1Si;O, melt. It is concluded that 


compositional variations and measurement techniques are 


causes for these differences. These problems have been 


reviewed recently by Day and Fenn (1982). 
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It 1s emphasized that the equimolar solubility of water 
is enhanced in increasingly peralkaline melts and this 
variation is not predicted by the model of Burnham (1979). 
The model of Burnham (1981) does reproduce the solubility 
determinations for PHO at 970 bars (wt% H,O calc. = 5. O02) 
and EA1.2K at 1620 bars (wt% H20 calc. = 4.68) but 
overestimates all the other melt solubilities in Table 8. 
Further discussion on solubilities in feldspathic melts is 
outside the scope of this study which pertains to 
haplogranitic melts and natural granitic melts. 

The fact that phonolitic melts show strongly enhanced 
solubility of H20 over S$i0O2-Saturated compositions indicates 
the potential for a far more dramatic role for H20 in the 
crystallization and emplacement of phonolitic magmas. 
Generally, it is clear that a range of water solubilities 
are to be expected in granitic melts of variable 
composition. These differences must be taken into account in 
the geochemical modelling of processes such as pegmatite 
formation, hydrothermal activity, crystallization, volatile 


exsolution, and eruptive mechanisms for granitic magmas. 


E. Geological application 
The experimentally determined solubility of water in 


haplogranitic melts and determinations of water in rhyolitic 


melt inclusions in phenocrysts may be used to estimate the 


minimum pressure of crystallization for the Bishop Tuff. 
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Druitt et a]. (1982) determined the concentration of water 
in rhyolitic glass in quartz from the basal Plinian airfall 
Soeeccitehe BishoprDuliftorbes4<9<+0E5% wAevacunmatcusion gas 
analyser method was used. The glass contained 77.6% Si0., 
13.3% Al203, 0.6% FeO, 0.5% MgO, 0.3% CaO, 3.9% Na-2O, and 
4.6% K20 (A. T. Anderson, personal communication). The gram 
molecular weight of the anhydrous melt (286 g) was 
calculated by the method of Burnham (1979). The Na/(Natk) 
and Si/Al ratios are 0.563 and 4.95, respectively, and near 
to or within the range investigated. The alkalies Na and K 
account for about 90 mole % of the exchangeable cations and 
the range of exchangeable cations to aluminum is 0.96. 
Therefore, the equimolar solubility for the stated 
composition is estimated to be X = 0.0111 x (£(H20),bars) . 
Accordingly, the pre-eruption fugacity of water estimated 
from the 4.9 +0.5% in the melt inclusions at 800°C is 1630 
+200 bars. The pressure of water is estimated to have been 
2100 +200 bars. The implied minimum depth of origin of these 
phenocrysts is 7.8 +1.1 km. The known presence in these 
inclusions of other components such as CO; (A. T. Anderson, 
personal communication) would increase the minimum pressures 
and depths of crystallization. Druitt et al. (1982) reported 


a minimum concentration of water (2 wt%) in rhyolitic 


inclusions in quartz phenocrysts from the last erupted Mono 


Lobe. Our solubility measurements suggest minimum fugacities 


and pressures of H20 of 480 and 530 bars, respectively. The 


direct measurements of pre-eruption water contents (Druitt 
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et al., 1982) and solubility of water in related 
haplogranitic melts at 1-2 kbar and 800°C (this study) yield 
more direct estimates of pressures than are available from 
mineral equilibria (Hildreth, 1979). 

The increase in water solubility with increasing 
peralkalinity observed in this work has at least one major 
implication for the petrogenesis of hydrous peralkaline 
melts. Silicic peralkaline volcanic suites often exhibit the 
chemical characteristic of increasing peralkalinity with 
increasing magmatic evolution. This effect was attributed by 
Bowen (1928) to the predominant crystallization of feldspar 
from such melts which results in residual peralkaline melts 
becoming more extremely peralkaline (i.e. the "plagioclase 
effect"). The data from this study show that one effect of 
increasing peralkalinity on an evolving magma system is to 
extend the percentage of crystallization which can occur 


before the residual magma becomes saturated with water. 
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VI. Major element partitioning 


A. Introduction 

Crystallization of hydrous granitic magmas at depth 
commonly results in the eventual saturation of the residual 
melt phase with water. Water saturation combined with 
further crystallization causes vesiculation of an aqueous 
fluid from the magma by the process of second boiling. The 
composition of the aqueous phase released in this manner is 
an meant parameter in alteration of the country rocks 
and in the reconstruction of the chemical compositions of 
magmas which yield holocrystalline granitic plutons. 

Chlorine-rich, magmato~hydrothermal systems are often 
host to economic concentrations of various metals (Burnham, 
1979) and, therefore, the composition of chlorine-bearing 
fluids in equilibrium with granitic melts has been 
extensively investigated by Burnham (1967 and data in Clark, 
1966). Pichavant (1981) has produced equivalent data for 
boron-bearing systems. Examples of fluorine-rich 
Magmato-hydrothermal systems occur in several localities 
(Bailey, 1977) illustrating the need for a better 
understanding of their geochemistry. This study was 


undertaken to characterize the chemical composition of the 


aqueous fluid phase in equilibrium with fluorine-rich 


granitic melts. 
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B. Experimental and analytical methods 

The starting compositions chosen for these experiments 
are shown in Figure 26. The glass compositions are centered 
on the 1 kbar haplogranite minimum melting composition 
(Tuttle and Bowen, 1958) and represent two compositional 
series ; KOR - EOR - NOR along the exchange vector KNa_,, 
and EAK1.5 - EAK1.2 - BOR - EALO.8 - EALO.5, along the 
exchange vector Al.2Na_,K_,. These compositions were chosen 
in order to evaluate the effect of melt Na/K ratio and 
alkali/aluminum ratio on melt/vapor partitioning. 

The glasses used in the hydrothermal runs were 
synthesized from reagent-grade Na2CO3, K2CO3, Al203 and SiO, 
by uSing standard techniques described by Bowen and Schairer 
(1956). The compositions of the anhydrous glasses used as 
Starting materials were determined by electron microprobe 
(Table 11). Experiments were conducted in cold-seal pressure 
vessels at 1 kbar and 800°C. Fifty milligrams of powdered 
glass was added to fifty milligrams of HF solution (1, 2 or 
4 wt% HF) in platinum capsules. The capsules were welded and 
checked for possible leakage by testing for weight loss 
after 15 minutes in an 110°C drying oven. Leak-free capsules 
were loaded into the cold-seal vessels and were raised to 
run pressure using an air-driven hydraulic pump. Pressurized 
capsules were heated to run temperature by using 
Kanthal-wound resistance furnaces mounted concentrically 
about the vessels along a vertical axis. Temperatures were 


monitored with Chromel-Alumel thermocouples and the actual 
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Fagure 26. Starting compositions in K20-Na20-Al,0;3;-Si0.-. 
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Table 11 
Glass Na20 
EOR 3.94 
NOR Aan 
KOR 2.41 
EALO.5 Zoe 
EALO.8 SUA sis 
BAK i a2 4.41 
EAK1.5 Suits’ 


K20 Al203 S102 


ey 11.48 80.0 
Sue S Wig6s pee 
6.3/7 hee 5 yy Orou 
2083 123.0 Salen 
Aoi Seis 78.4 
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-uncertainties quoted as percent of the amount present 
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: Microprobe analyses of starting glasses. 
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temperatures are believed accurate to +5°C. Pressures were 
measured with Bourdon-tube gauges previously calibrated with 
a Heise gauge and are accurate to + 20 bars. Run durations 
were 10° seconds for all runs except those which form a time 
series test of equilibrium (runs 12b, 16, 27 and 35). The 
charges were quenched isobarically by opening the vessel to 
the pressure line and by pumping for the duration of the 
quench . 

Reversal experiments were accomplished using the 
fluorine-bearing glasses from forwards experiments. 
Fluorinated glasses were ground in an agate mortar, rinsed 
with distilled water, and fifty milligrams of glass were 
added to fifty milligrams of triple-distilled water. 
Reversals were then run under equivalent conditions as the 
forwards experiments. As a check on the inert behavior of 
platinum in the presence of fluorine-bearing aqueous fluids 
at this pressure and temperature, duplicate experiments were 
conducted using gold capsules and equivalent results were 
obtained (Fig. 27). 

Run capsules which suffered no weight change (< 1 
milligram) and no visible signs of breach during the run, 
were ultrasonically cleaned in acetone and frozen in liquid 
nitrogen prior to opening. The frozen capsules were slit 
longitudinally with a razor blade, immersed in 
triple-distilled 2N HNO;, and returned to the ultrasonic 
bath to help dissolve quench fluorides. Atomic absorption 


Spectrometric analyses were carried out using the Dept. of 
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mole % metals 


Run# Details 
12b 250,000 sec (5% days) 
16 10,000,000 sec (115 days) 


AU capsule 
35 2,500,000 sec (28% days) 
1,000,000 sec reversal 
1,000,000 sec reversal 
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Geology Perkin Elmer 503 Spectrophotometer. Standards for 
Na, K, Al and Si were prepared in HNO; media. For certain 
experiments (runs 16, 20, 27 and 29) an aliquot of the 
remaining fluid was titrated to neutral PH by using 8% NaOH, 
combined with a total ionic strength buffer, and analyzed 


for fluorine by using an Orion fluorine specific electrode. 


C. Results 

The results of Na, K, Al and Si analyses of the fluid 
phase are presented in Table 12. Runs # 12b, 16, 27 and 35 
form a time series test for equilibrium and runs # 54 and 55 
are reversals. Run # 35 is a gold-encapsulated experiment. 
The attainment of chemical equilibrium in the partitioning 
of Na, K, Al and Si between vapor and melt is evident in the 
results of Table 12. Figure 27 is a projection of the solute 
contents, of the vapor phase into molar Na-K-Al-Si space. 
Figure 27 indicates the reproducibility obtained in the 
analysis, of thexsolute; andethese.projeetions are used in 
subsequent figures. The largest scatter in the data is in Al 
and Si content due to the relatively poor precision of the 
atomic absorption analyses at these low levels (actual 
analyzed solutions contained 1-10 ppm Al and 10-50 ppm Si). 
The solute yields for EOR runs with 1, 2, and 4wt% fluorine 
and chlorine added as HF and HCl, respectively, are shown in 
Figure 28. Fluorine and chlorine have similar effects on the 
Solubility of Na and K. Na and K contents increase rapidly 


with added F or Cl. The dissolved contents of Si, and 
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Figure 28. Solute yields for Na, K, Al and Si. 
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PapiLevsi2 +. Na, Ka’) ana: <i analyses of the vapor phase. 


Run Composition 


2) EOR2F 
1 2DtS BOR2E 
16 EOR2F 
2g EOR2F* 
55 EOR 2F 


54 |. BOR2F*x 
55 EOR2F xx 


SZ NOR2F 
33 KOR2F 
20 EOR4F 
23 EORIF 


50 BEAK 1.2-2F 
Sil EALO.8-2F 
72 EALO.5-2F 
74 BAK. Once 


15 EOR2C1 
22 EOR1C1l 
9 EOR4C1 


1 
2 
10 
1 
25 
1 


] 
1 
1 
1 


1 
; 
1 
1 
1 


1 
] 
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Duration Na 
(sec) 


,000,000 425 
POUU ROU 450 
,000,000 443 
,000,000 366 
,000,000 490 
,000,000 393 
,000,000 RE 
,000,000 437 
,000,000 215 
, COO O00"* M966 
,000,000 Zn9 
,000,000 835 
,000,000 528 
,000,000 402 
O00, 000. W071 


,000,000 1940 
,000,000 250 
20007000 3510 


-* gold capsule experiment 


-** reversal experiment 
- all analyses quoted as ppm; 
as follows: Na and K (+25 ppm); Al 
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especially Al, increase Substantially with increasing 
fluorine in the vapor phase, while chlorine has little or no 
effect on Si solubility and an apparent negative effect on 
Piesolubility, 

Yields of fluorine from the vapor phase are presented 
in Figure 29. A vapor/melt partition coefficient of 0.08 
holds for the range 2 to 4 wt% F. 

Figure 30 presents the molar projection of melt and 
vapor compositions for these experiments. Tie lines join 
melt compositions (Solid triangles) to coexisting fluid 
compositions (open triangles). Clearly, the vapor 
composition is strongly dependent on melt composition. These 
projections show that all vapor compositions are more sodic 
and alkalic than the corresponding melt compositions. The 
alkali/aluminum ratios of the solute are directly controlled 
by the melt composition. 


For convenience, vapor/melt distribution coefficients 


may be defined as follows 


D(Na/K) = (Xya/Xk) vapor and D(AK/Al) = (Xyaak/¥ar) vapor 
(Xy/Xq) melt (Xnatk/Xa1) melt 


The values of these distribution coefficients for the melts 
studied are given in Table 13. The value of D(Na/K) ranges 
from 1.9 to 2.4 and decreases with increasing peralkalinity, 
while the value of D(Ak/Al) ranges from 2.3 to 5.4 and is 


highest for peralkaline melts. 
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Figure 29. Solute yields for fluorine. 
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Figure 30. Experimental results of the partitioning 
experiments; effect of melt composition. 
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Table 13 : Melt/vapor distribution coefficients 


Melt Composition D (Na/K) D (Ak/Al) 
EALO.5 2.4 Sil) 
EALO.8 2.4 20 
EOR 20 De0 
EAK1.2 ARS, 4.7 
EAK1.5 ee. 5.4 
NOR 2.4 Sl) 
KOR 2.0 25% 


-uncertainties in melt/vapor partition coefficients are; 


D(Na/kK) + 40% ; D(Ak/Al) 270%. 


(fANQAA) 


i 


Nu 


ss + 
; 


met 


_ 
a 


a i, ae 


fi 


D. Discussion 


Although the data presented above are not extensive 
enough to permit modelling of the exact stoichiometry of 
fluoride complexation in these fluids, the following 
observations are significant. 

The solubility of Al is extremely sensitive to fluorine 
concentration, implying the existence of some form of 
aluminofluoride species in solution at this pressure and 
temperature. It cannot be concluded, however, that the Al 
solubility increase is due entirely to the formation of a 
Simple Al-F complex. Increased Na, K and Si solubility with 
increasing fluorine concentration may be indicating the 
Stabilization of one or more alkali-aluminofluoride 
complexes (cryolite-like species?) in solution. What does 
seem certain, however, is that fluorine does not compete 
with Al for a large proportion of the dissolved alkalies to 
the same extent as Cl. Anderson and Burnham (1983) have 
proposed that the low solubility of Al in hydrothermal 
fluids containing Cl is a direct consequence of alkali 
chloride complexation which decreases the amount of 
dissolved alkalies available for complexation in 
alkali-aluminosilicate complexes (feldspar stoichiometry). 
Decreasing Al solubility as a function of chlorine 
concentration in Figure 28 probably results from this 
effect. 

Another aspect of fluoride complexation is that the 


effects of fluorine on hydrothermal fluids is remarkably 
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large for relatively low concentrations of dissolved Pac cAt 
equivalent molalities of HF and HCl in these experiments, 
the vapor phase in chlorine-bearing experiments contains an 
order of magnitude more dissolved halide than the 
corresponding fluorine-bearing experiment, due to the 
different vapor/melt partition coefficients for fluorine 
(Hards, 1978) and chlorine (Kilinc and Burnham, 1972) 
Comparison with the data of Burnham (1967) for the 
Spruce Pine pegmatite melt-HC1-H20 is difficult due to the 
higher pressure range of Burnham's study (2-6 kbars); 
however, the pressure dependence of the composition of the 
chlorine-bearing fluid is not large (Fig. 2.7 in Burnham, 
1967). Figure 31 compares the compositions of 
fluorine-bearing fluids coexisting with haplogranitic melts 
(this study) with the results of Burnham (1967) for 
chlorine-bearing fluids and Pichavant (1981) for 
boron-bearing fluids. From this comparison it is concluded 
that the solute contents of F-bearing fluids exhibit higher 
Al/(NatK) ratios than either Cl- or B-bearing fluids while 


Cl-bearing fluids exhibit lower Si/(Na+K) ratios than either 


F- or B-bearing fluids. 


E. Geological implications 


The physical and chemical conditions which permit the 


coexistence of granitic magmas and aqueous fluids have 


received extensive consideration in the literature (Burnham, 


1967, 1979; Holland, 1972; Burnham and Ohmoto, 1980). 
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Figure 31. Comparison of the compositions of 
fluorine-bearing, chlorine-bearing and boron-bearing 
fluids coexisting with granitic melts (see text for data 
sources). 
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An aqueous fluid phase may be evolved from a granitic 
magma in two ways; decompression and crystallization. The 
rapid ascent of a water-undersaturated melt is accompanied 
by a pressure decrease which, in turn, decreases the 
solubility of water in the melt and results in vesiculation. 
Alternatively, the crystallization, at constant total 
pressure, of dominantly anhydrous and water-poor phases from 
such a melt will also result in the water saturation of the 
diminishing melt volume due to the increasing proportion of 
dissolved water remaining in the melt (i.e., second, 
resurgent, or retrograde boiling). 

One or both of these processes are involved in every 
geological situation where magmatic water forms as a 
separate phase. If the separation of an aqueous phase from 
an ascending magma occurs at depths which correspond to 
pressures greater than the critical region of the aqueous 
fluid then the immediate result of aqueous phase separation 
is a supercritical fluid. The present study simulates the 
case of a supercritical aqueous fluid which is equilibrated 
with such a melt. 

If the aqueous phase separates at shallow, sub-volcanic 
depths corresponding to pressures of several hundred bars 
then the H,O-NaCl critical curve may intersect the solidus. 
In this case the aqueous phase may separate as two distinct 
fluid phases, a vapor and a brine, ise. critical behavior: 
the coexistence of sub-critical brine and melt 


Evidence for 


is provided by Roedder and Coombs (1967). Primary boiling 
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also occurs if magmatic aqueous fluids produced at depth 
ascend into the P-T region of critical behavior, (or if high 
temperature, low pressure fluids from more mafic, 
subvolcanic melts cool into this P-T region. 

The formation of deep-seated, intrusion-hosted 
pegmatites is an example of the separation of a 
Supercritical fluid phase and subsequent precipitation of 
dissolved species in the absence of boiling (Jahns and 
Burnham, 1969) Many pegmatitic greisens may form in this 
manner (Shcherba, 1970, Burt, 1981). In contrast, the fluid 
inclusion evidence of simple or first boiling is widespread 
in shallow, subvolcanic paemmeoeh Geet necee systems such as 
porphyry Cu deposits (Burnham and Ohmoto, 1980; Gustafson, 
1978; Titley and Beane, 1980). 

Fluorine has a high melt/fluid partition coefficient 
compared with chlorine ( Hards, 1978; Kilinc and Burnham, 
1972) . Therefore, hydrous granitic magmas, containing 
similar amounts of fluorine and chlorine, which become 
water-saturated will yield Cl-dominated aqueous fluids. 
Fluid inclusions often record fluid salinity and 
compositions which reflect post-solidus modification by 
primary boiling. 

Despite the affinity of chlorine for aqueous Eluids, 
evidence for mild to extreme enrichments of fluorine in 
certain granitic magmas and associated fluids is evidenced 


by petrochemical Studies of silicic hypabyssal and volcanic 


rocks (Kovalenko, 1973; Bailey, 1977: Eadington and Nashar, 
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1978; Burt et al., 1982; Christiansen et al., 1983, 1984) 
and the mineral chemistry of equivalent plutons (Pauly, 
1960; Gunow et a]., 1980: Nedachi, 1980; Czamanske eft a/., 
1981). Fluorine enrichment in such Systems as deep-seated 
pegmatitic greisens (Burt, 1981), shallow Mo-porphyries, 
stockwork Mo and vein-type Sn and W deposits (Westra and 
Keith, 1981; Huspeni et a]., 1984) and lithophysal topaz 
Ehyolttes (Burt Ef%alv7 1982) Christiansen sere al-+) + 1983+ 
1984) all clearly indicate that F-dominated 
magmato-hydrothermal systems can and do occur. The causes of 
such fluorine concentration are not clearly understood and 
several factors may contribute to their formation. The most 
Crucial factor may be the early depletion of the systems in 
chlorine by fluid escape. Certainly, the evidence for 
multiple intrusion and protracted magmatic evolution of 
fluorine-rich systems does exist (White ef a]., 1980). 
Alternatively, the magmas may have originated in Cl-poor 
source regions due to prior metasomatic and/or magmatic 
processes which resulted in chlorine depletion. 

Regardless of the causes for the development of 
fluorine-rich systems, one of the most important 
considerations for the magmato-hydrothermal systems which 
they yield is the role of fluorine in enhancing the 
solubilities of dissolved metals. In most cases, fluorine is 
eventually partitioned into a fluid phase which is 


responsible for its mobility and concentration. 
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F-rich systems are also characterized by evidence of Al 
mobility, a feature absent in Cl-dominated systems (Titley 
and Beane, 1980: Anderson and Burnham, 1983; Rubie and 
Gunter, 1983). The results of this Study provide an 
experimental confirmation of the mobility of Al in close. 
association with fluorine-rich igneous rocks as evidenced by 
distal fluor-muscovite veins in fluorine-rich skarns 
(Barton, 1982); massive cryolite roof rocks (Pauly, 1960); 
topaz-lined vesicles in lithophysal rhyolites (Burt et al., 
1982; Christiansen et a/]., 1983, 1984), and cryolite + 
Fluorite + elpasolite daughter crystals in fluid inclusion 
of hypabyssal ongonites (Naumov et a/]., 1971, 1977). 

Unfortunately, the results presented here can only be 
applied qualitatively at this time because the melt-fluid 
system investigated is multivariant in the absence of 
calibrated buffer phases for all components except hydrogen. 
Quantification of fluorine speciation and fluoride complex 
Stability in aqueous fluids may be measured by mineral 
solubility studies (Anderson and Burnham, 1967; Frantz et 
al., 1981) and exchange mineral/fluid exchange equilibria 
(Orville, 1963; Barton and Frantz, 1983). The thermodynamic 
data from these experiments will eventually allow 
quantitative modelling of aqueous fluid behavior . Extension 
of these studies to mineral + melt + fluid systems should 
eventually permit further thermodynamic modelling of 


water-saturated natural Silicate melts. 
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VII. Summary and conclusions 


A. Summary of results 

Various aspects of the solution of fluorine in 
aluminosilicate melts have been investigated in this thesis. 
The following is a short summary of the results. 

The viscosity of polymerized melts in Naz0-Al,03-Si0, 
is strongly reduced by the addition of fluorine. These 
viscosity reductions are similar to the effect of added 
water on melt viscosity. The viscosity-reducing effect of 
fluorine increases with SiO, in melts along the join 
NaAlO2-SiO2 and at 75 mole percent S$i02, fluorine-bearing 
melts increase in viscosity from peralkaline through albite 
to peraluminous compositions. Due to the lower activation 
energy of viscous flow in fluorine-bearing melts, the 
viscosity reductions with fluorine addition increase with 
decreasing temperature. 

The chemical diffusion of fluorine in these melts is 
much faster than the self-diffusion of oxygen but 
significantly slower than the chemical diffusion of water. 
The chemical diffusion of .fluorine is a binary, 
concentration-independent, composition-dependent exchange 
with oxygen. This behavior contrasts with the chemical 
diffusion of water which may involve alkali exchange. 
Volatilization of fluorine from dry aluminosilicate melts is 


predominantly as SiF, and this property has been used to 
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investigate fluorine diffusivity at 1 atmosphere. Fluorine 
increases the diffusivity of oxygen and network-modifying 
cations in aluminosilicate melts as does water. The pressure 
dependence of F-O interdiffusion is very Slight in the 
pressure range of 10 to 15 kilobars but from 0.001 to 10 
kilobars an increase in interdiffusivity is observed. F-O 
interdiffusivity decreases along the join NaAl0O,-Si0O, and, 
at 75 mole percent SiO2, F-O interdiffusion is slower in 
albite melt than in peralkaline and peraluminous melts. 

The solubility of water in haplogranitic melts at 1 to 
2 kilobars and 800°C is greater in peralkaline and 
peraluminous melts than in the 1 kilobar granitic minimum 
melt. Phonolitic melts dissolve twice as much water as 
gGranitic melts. Fluorine slightly decreases the solubility 
of water in granitic melts. 

Fluorine, in contrast with chlorine, increases the 
solubility of aluminum in magmato-hydrothermal fluids. The 
distribution of Na and K between granitic melts and 
coexisting aqueous fluids is controlled by the melt Na/K 
ratio and is similar for fluorine- and chlorine-bearing 
aqueous fluids. Despite the strong preference of fluorine 
for the melt phase, fluorine increases the aqueous 
solubility of each of Na, K, Al and Si in the fluid phase. 


The alkali/aluminum ratio of the fluid phase is influenced 


by the melt alkali/aluminum ratio ; however, the 


alkali/aluminum (molar) ratio of the fluorine-bearing 


aqueous fluid phase is greater than 1 in all experiments. 
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B. Conclusions 


Melt structure 

The dissolution of fluorine in dry aluminosilicate 
melts involves the replacement of oxygen in SP=O=(A DIST) 
bridges by Si-F or Al-F units. Thus, fluorine depolymerizes 
the three-dimensional structure of polymerized melts in 
Na20-Al203;-Si0O2. The viscous flow of aluminosilicate melts 
with added fluorine or water is similar, indicating that 
(Al,Si)-F and (Al,Si)-OH units within the melt have similar 
resistance to viscous flow. Despite this similarity, the 
transport of fluorine within aluminosilicate melts is unlike 
that -of water. Fluorine transport occurs by an effective 
binary exchange with oxygen, whereas water diffusion 
probably involves alkali exchange. Failure of the 
Stokes-Einstein relationship to predict | 
viscosity-diffusivity relationships in fluorine-bearing 
melts suggests that the mechanisms and/or species involved 
in viscous and diffusive transport in these melts are not 
equivalent. 

In hydrous aluminosilicate melts, fluorine may compete 
directly with water for dissolution sites, thus decreasing 
the solubility of water in these melts. Fluorine is very 
effective in decreasing the activity of feldspar-forming 
components in hydrous aluminosilicate melts (Manning, 1981) 


and the decrease in water solubility with fluorine addition 
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may be an indirect consequence of the decrease in the 
activity of NaAlO, units. The solubility of water in 
fluorine-free melts is strongly influenced by the 
alkali/aluminum ratio but not by the Na/K ratio. The 
behavior of water solubility as a function of 
alkali/aluminum ratio strongly suggests the operation of 
more than one mechanism of water solution. The preference of 
fluorine for Al-bearing species in supercritical aqueous 
fluids indicates that fluorine does not compete strongly 
with aluminum for the coordination of dissolved alkalies as 


is the case for chlorine. 


Geological implications 

Silicic volcanic rocks which contain up to 2 weight 
percent fluorine are often characterized by low eruption 
temperatures, extensive lava flows, crystal-poor glasses, 
and topaz-lined lithophysae (Burt et a]., 1982). Such 
volcanics often have trace element signatures consistent 
with an origin involving extreme degrees of crystal 
fractionation (Christiansen et a]., 1984). The subvolcanic 
granites of equivalent chemistry are often host to economic 
concentrations of several elements whose enrichment can be 
explained by accumulation in residual melts during 
protracted magma fractionation (Christiansen et a]., 1983). 

Several of these characteristics of fluorine-rich 
s are readily explicable on the basis of the 


igneous rock 


results of this and previous studies on the role of fluorine 


Je a 
(688! ,. (a 06 nepngh gaat 


>| 


su S62! au wiresacs 
(jay 1 eo IEG Choi 
' ah 
£ ee: dog* kd sete 3 


(se 4 SHe! 


eetiat: 
fox s~strtnoutt i< 298m 
4 ais 7 
| o,?. 


a6s So slasa aft os) 


ts 


in aluminosilicate melts. 

Perhaps the most important consequence of fluorine 
concentration in late stage igneous melts is that fluorine 
can appreciably extend the magmatic "lifespan" of 
crystallizing granitic magmas by depressing the 
liquidus/solidus phase relationships several hundred 
degrees. This depression occurs in both dry (Koster van 
Groos and Wyllie, 1968) and wet (Manning, 1981) systems and 
one can reasonably infer, therefore, that the effect extends 
over the range of water-undersaturated conditions pertinent 
to the formation of many fluorine-rich silicic volcanics. 
Fluorine has a high melt/vapor partition coefficient (Hards, 
1978) and appears to behave as an incompatible element in 
many silicic volcanic systems (Bailey, 1977). Progressive 
crystallization in such systems leads, therefore, to 
fluorine concentration in residual melts which in turn 
experience decreasing solidus temperatures. The fact that 
water-undersaturated, fluorine-rich rhyolites erupt at 
temperatures in the range 630 to 850°C (Bikun, 1980) is 
clear evidence of the influence of fluorine in such systems. 

The influence of fluorine on melt diffusivities 
probably increases the nucleation and growth rates of 
observed phenocryst phases in fluorine-rich rhyolites. The 
effect is the "mineralizing action" that Buerger (1948) 
attributed to both F~ and OH. The 
concentration-independence of fluorine diffusion means that 


the effect of fluorine on melt diffusivities should be felt 
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at very low concentrations of dissolved fluorine, and that 
local enrichments of fluorine in the vicinity of growing 
crystals is unlikely. If protracted crystal fractionation is 
Significant in the petrogenesis of fluorine-rich rhyolites, 
then crystal fractionation must operate as an efficient 
process in these melts. The viscosity-reducing effect of 
fluorine may be an important parameter that separates the 
crystal fractionation behavior of fluorine-rich rhyolites 
from that of fluorine-free rhyolites. Lower viscosities also 
result in the formation of large individual lava flows which 
travel considerable distances from their volcanic vents. 
Finally, it 1s clear that some crystal-poor 
fluorine-rich magmas are sampled volcanically (Burt et al., 
1982). The result is rhyolites whose vesicles are lined with 
the quench products of an aqueous fluid phase which was 
released during eruption. Quench alkali fluorides in the 
lithophysae are rapidly leached by groundwaters, but the 
relatively insoluble fluor-topaz remains as evidence of the 
aluminum-carrying capacity of fluorine-bearing aqueous 


fluids. 
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